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Abstract 

This paper outlines astrophysical issues related to the long term fate of the universe. We consider the 
evolution of planets, stars, stellar populations, galaxies, and the universe itself over time scales which greatly 
exceed the current age of the universe. Our discussion starts with new stellar evolution calculations which 
follow the future evolution of the low mass (M type) stars that dominate the stellar mass function. We 
derive scaling relations which describe how the range of stellar masses and lifetimes depend on forthcoming 
increases in metallicity. We then proceed to determine the ultimate mass distribution of stellar remnants, 
i.e., the neutron stars, white dwarfs, and brown dwarfs remaining at the end of stellar evolution; this 
aggregate of remnants defines the “final stellar mass function”. At times exceeding ^1-10 trillion years, 
the supply of interstellar gas will be exhausted, yet star formation will continue at a highly attenuated 
level via collisions between brown dwarfs. This process tails off as the galaxy gradually depletes its stars by 
ejecting the majority, and driving a minority toward eventual accretion onto massive black holes. As the 
galaxy disperses, stellar remnants provide a mechanism for converting the halo dark matter into radiative 
energy. Posited weakly interacting massive particles are accreted by white dwarfs, where they subsequently 
annihilate with each other. Thermalization of the decay products keeps the old white dwarfs much warmer 
than they would otherwise be. After accounting for the destruction of the galaxy, we consider the fate of 
the expelled degenerate objects (planets, white dwarfs, and neutron stars) within the explicit assumption 
that proton decay is a viable process. The evolution and eventual sublimation of these objects is dictated 
by the decay of their constituent nucleons, and this evolutionary scenario is developed in some detail. After 
white dwarfs and neutron stars have disappeared, galactic black holes slowly lose their mass as they emit 
Hawking radiation. This review finishes with an evaluation of cosmological issues that arise in connection 
with the long-term evolution of the universe. We devote special attention to the relation between future 
density fluctuations and the prospects for continued large-scale expansion. We compute the evolution of 
the background radiation fields of the universe. After several trillion years, the current cosmic microwave 
background will have redshifted into insignificance; the dominant contribution to the radiation background 
will arise from other sources, including stars, dark matter annihilation, proton decay, and black holes. 
Finally, we consider the dramatic possible effects of a non-zero vacuum energy density. 
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I. INTRODUCTION 


The long term future of the universe and its contents is a topic of profound scientific and philosophical 
importance. With our current understanding of physics and astrophysics, many of the questions regarding 
the ultimate fate of the universe can now be quantitatively addressed. Our goal is to summarize and 
continue the development of a quantitative theory of the future. 

Investigations of the early universe at both accessible and inaccessible energies have become com¬ 
monplace, and a great deal of progress within this discipline has been made (see, e.g., Weinberg, 1972, 
1977; Kolb & Turner, 1990; Linde, 1990; Peebles, 1993; Zuckerman & Malkan, 1996). On the other hand, 
relatively little work has focused on the future of the universe. The details of the fiery denouement in 
store for a closed universe have been outlined by Rees (1969), whereas an overview of the seemingly more 
likely scenario in which the universe is either open or flat, and hence expands forever, was set forth in the 
seminal paper Time Without End (Dyson, 1979). The development of an open universe was also considered 
in detail by Islam (1977, 1979). The spirit of Rees, Islam, and Dyson’s work inspired several follow-up 
studies (see also Rees, 1981). The forthcoming evolution of very low mass stars has been discussed in gen¬ 
eral terms by Salpeter (1982). The effects of matter annihilation in the late universe were studied (Page 
& McKee, 1981ab), and some aspects of proton decay have been explored (Dicus et ah, 1982; Turner, 
1983). Finally, the possibility of self-reproducing inflationary domains has been proposed (Linde, 1988). 
In general, however, the future of the universe has not been extensively probed with rigorous calculations. 

Because the future of the universe holds a great deal of intrinsic interest, a number of recent popular 
books have addressed the subject (e.g., Davies, 1994; Dyson, 1988; Barrow & Tipler, 1986; Poundstone, 
1985). Authors have also grappled with the intriguing prospects for continued life, both human and 
otherwise, in far future (e.g., Dyson, 1979; Frautschi, 1982; Barrow & Tipler, 1986; Linde, 1988, 1989; 
Tipler, 1992; Gott, 1993; Ellis & Coule, 1994). Our aim, however, is to proceed in as quantitative a 
manner as possible. We apply known physical principles to investigate the future of the universe on 
planetary, stellar, galactic, and cosmic scales. The issue of life, however alluring, is not considered here. 

In standard Big Bang Cosmology, evolutionary epochs are usually expressed in terms of the redshift. 
When considering the far future, however, time itself is often the more relevant evolutionary measure. The 
immense dynamic range of time scales r involved in the subject suggests a convenient logarithmic unit of 
time rj, defined by 


V = logio 


L(lyr)J • 


( 1 . 1 ) 


We refer to a particular integer value of ry as a “cosmological decade”. For example, the current age of the 
universe corresponds to ry « 10. 

The article of faith inherent in our discussion is that the laws of physics are constant in time, at least 
over the range of time scales 10 < ry < 100 under consideration. There is no general guarantee that this 
assumption holds. Nevertheless, modern cosmology suggests that physical laws have held constant from 
the Planck time to the present, i.e., over cosmological decades spanning the range —50 < ry < 10, and there 
is little reason to expect that they will not continue to do so. We also implicitly assume that all of the 
relevant physics is known (with full awareness of the fact that our version of the future will be subject to 
revision as physical understanding improves). 

This paper is organized in roughly chronological order, moving from events in the relatively near future 
to events in the far future. In section §11, we discuss physical processes that affect conventional stellar 
evolution; these processes will take place in the time range 10 < ry < 15. In §III, we discuss events which 
lead to the disruption and death of the galaxy; these processes unfold over a time range 15 < ry < 25. 
Marching further into time, in §IV, we discuss the fate of stellar objects in the face of very long term 
processes, including proton decay (30 < ry < 40), and Hawking radiation (60 < ry < 100). In §V, we 
broaden our scope and focus on the long term evolution of the universe as a whole. We conclude, in §VI, 
with a general overview of our results. Since physical eschatology remains embryonic, we emphasize the 
major unresolved issues and point out possible avenues for further research. 
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II. THE END OF CONVENTIONAL STELLAR EVOLUTION 


At the present epoch, stars are the cornerstone of astrophysics. Stars mediate the appearance and 
evolution of galaxies, stars are responsible for evolving the chemical composition of matter, and stars 
provide us with much of the information we have regarding the current state of the universe. 

For the next several thousand Hubble times, conventionally evolving stars will continue to play the 
central role. We thus consider the forthcoming aspects of our current epoch, which we term the Stelliferous 
Era. In particular, the fact that the majority of stars have barely begun to evolve motivates an extension of 
standard stellar evolution calculations of very low mass stars to time scales much longer than the current age 
of the universe. We also discuss continued star formation within the galaxy, and the final mass distribution 
of stellar remnants. 


A. Lifetimes of Main Sequence Stars 


Low mass stars are by far the most commonplace (e.g., Henry, Kirkpatrick, & Simons, 1994), and they 
live for a long time. To a working approximation, the main sequence (core hydrogen burning) lifetime of 
a star depends on its mass through the relation 


T* 


10^°yr 


r M* 1 

-IMq. 


— a 


(2.1a) 


where the index a « 3 — 4 for stars of low mass. In terms of cosmological decades ry, we obtain 


ry* = 10 - a log^o [AT*/IM©]. 


(2.1b) 


Thus, for example, ry* « 13 for a small star with M* = 0.1 M©. Indeed at the present time, only stars 
with masses M* > 0.8M© have had time to experience significant post-main sequence evolution. Hence, a 
large fraction, 


/ = 


(dN/dm)dm 

J IVlmin ^ _ ' _ _ _ 

j^2Z^‘^N/dm)dm 


80%, 


( 2 . 2 ) 


of all stars ever formed have yet to experience any significant evolution (here, dN/dm is the mass distri¬ 
bution - see §11.E). We are effectively still in the midst of the transient initial phases of the stelliferous 
epoch. 


Very little consideration has been given to the post-main sequence development of stars which are 
small enough to outlive the current age of the universe. An essay by Salpeter (1982) contains a qualitative 
discussion regarding the evolution of M stars (especially with respect to ^He production) but detailed 
stellar evolutionary sequences have not been presented in the literature. Nevertheless, there is a sizable 
collection of papers which discuss the pre-main sequence and main sequence properties of very low mass 
stars (e.g., Kumar, 1963; Copeland, Jensen & Jorgensen, 1970; Grossman & Graboske, 1971; D’Antona & 
Mazzitelli, 1985; Dorman, Nelson & Chau, 1989). The best comprehensive family of models spanning the 
M dwarfs and brown dwarfs is probably that of Burrows et al. (1993). Those authors devote attention 
to the formative cooling phases, as well as the exact mass of the minimum mass star (which for their 
input physics occurs at M* = 0.0767 M©). Evolution beyond 20 billion years was not considered (see also 
Burrows & Liebert, 1993). 

The dearth of information regarding the fate of the M dwarfs has recently been addressed (Laughlin, 
Bodenheimer, & Adams, 1996). We have performed a detailed series of stellar evolution calculations which 
follow the pre main-sequence through post main-sequence evolution of late M-dwarfs, yielding the following 
picture of what lies in store for the low mass stars. 

Newly formed stars containing less mass than M* ~ 0.25M© are fully convective throughout the 
bulk of their structure. The capacity of these stars to entirely mix their contents has several important 
consequences. First, these late M stars maintain access to their entire initial reserve of hydrogen, greatly 
extending their lifetimes in comparison to heavier stars like the sun which see their fuel supply constricted by 
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stratified radiative cores. Second, as recognized by Salpeter (1982), full convection precludes the buildup 
of composition gradients which are ultimately responsible (in part) for a star’s ascent up the red giant 
branch. The lowest mass stars burn all their hydrogen into helium over an ry = 13 time scale, and then 
quietly fade from prominence as helium white dwarfs. This general evolutionary scenario is detailed in 
Figure 1 (adapted from Laughlin et ah, 1996), which charts the path in the Hertzsprung-Russell diagram 
followed by low mass stars of several different masses in the range O.OSMq < M* < O.25M0. 

Upon emerging from its parent cloud core, the lowest mass star capable of burning hydrogen (M* « 
O.OSMq) descends the convective Hayashi track and arrives on the main sequence with a luminosity L* ~ 
IO-^Lq. The main sequence phase is characterized by gradual prolonged increase in both luminosity and 
effective surface temperature as hydrogen is consumed. Due to the relatively low prevailing temperature 
in the stellar core (Tc « 4 x 10® K), the proton-proton nuclear reaction chain is decoupled from statistical 
equilibrium, and the concentration of ®He increases steadily until 77=12.6, at which time a maximum mass 
fraction of 16% ®He has been attained. As the initial supply of hydrogen is depleted, the star heats up 
and contracts, burns the ®He, increases in luminosity by a factor of 10, and more than doubles its effective 
temperature. After ~11 trillion years, when the star has become 90% ^He by mass, a radiative core finally 
develops. The evolutionary time scale begins to accelerate, and hydrogen is exhausted relatively quickly in 
the center of the star. When nuclear burning within the modest resulting shell source can no longer provide 
the star’s mounting energy requirements, the star begins to contract and cool and eventually becomes a 
helium white dwarf. Stars with masses up to ~0.20 Mq follow essentially this same evolutionary scenario. 
As stellar mass increases, radiative cores develop sooner, and the stars perform increasingly dramatic 
blueward excursions in the H-R diagram. 

A star with a slightly larger mass, M* = 0.23 M©, experiences the onset of a radiative core when the 
hydrogen mass fraction dips below 50%. The composition gradients which ensue are sufficient to briefly 
drive the star to lower effective temperature as the luminosity increases. In this sense, stars with mass 
M* = 0.23Mq represent the lowest mass objects that can become conventional “Red Giants”. At these 
low masses, however, the full giant phase is not completed. Stars with initial mass M* < 0.5 Mq will be 
unable to generate the high central temperatures (Tc ~ 10® K) required for the helium flash; these stars 
abort their ascent up the giant branch by veering to the left in the H-R diagram in the manner suggested 
by Figure 1. 

The steady luminosity increases experienced by aging M dwarfs will have a considerable effect on the 
mass to light ratio of the galaxy. For example, as a 0.2 Mq star evolves, there is a relatively fleeting epoch 
(at 77 « 12) during which the star has approximately the same radius and luminosity as the present day 
sun. Given that M dwarfs constitute the major fraction of all stars, the total luminosity of the galaxy will 
remain respectably large, Tgai ~ IO^^Tq, at this future date. This luminosity is roughly comparable to the 
characteristic luminosity C* = 3.4 x IQio Lq displayed by present day galaxies (Mihalas & Binney, 1981). 


B. Forthcoming Metallicity Effects 

The foregoing evolutionary calculations assumed a solar abundance set. In the future, the metallicity 
of the galaxy will steadily increase as stars continue to process hydrogen and helium into heavy elements. 
It is thus useful to determine the effects of these metallicity increases. 


1. Stellar Lifetimes vs Metallicity 

First, it is possible to construct a simple scaling relation that clarifies how stellar lifetimes r* depend 
on the metallicity Z. The stellar lifetime is roughly given by amount of fuel available divided by the rate 
of fuel consumption, i.e., 

r* ~ M* A /T , (2.3) 

where M* is the stellar mass and X is the hydrogen mass fraction. For relatively low mass stars, the 
luminosity L obeys the scaling relation 

L ~ , (2.4) 
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where /r is the mean molecular weight of the star and where kq is the constant of proportionality appearing 
in the usual opacity relation for stars (Clayton, 1983). Thus, for a given stellar mass M*, the lifetime scales 
according to 


Ti, ^ KqX ^ 


(2.5) 


To evaluate the stellar lifetime scaling relation, one needs to know how the parameters kq, X, and fj, 
vary with metallicity. The opacity constant kq is roughly linearly dependent on the metallicity, i.e., 

kq ^ Z . (2.6) 

The mean molecular weight /i can be approximately written in the form 


{l + iX + Y/2) ’ ^ ^ 

where Y is the helium mass fraction (e.g., see Clayton, 1983). By definition, the mass fractions obey the 
relation 

X + Y + Z = 1. (2.8) 

Finally, for this simple model, we write the helium abundance Y in the form 

Y = Yp + fZ, (2.9) 


where Yp is the primordial abundance and the factor / accounts for the increase in helium abundance as 
the metallicity increases. Big Bang nucleosynthesis considerations indicate that Yp « 1/4 (Kolb & Turner, 
1990), whereas / ~ 2 based on the solar enrichment in Y and Z relative to the primordial values. Combining 
the above results, we obtain a scaling relation for the dependence of stellar lifetimes on metallicity, 

Th. ~ Y(1 - aZ) (1 - 6Y)^(2.10) 

where we have defined constants a = 4(l + /)/3 « 4 and b = 8/9 + 20//27 « 64/27. This result implies that 
stellar lifetimes have a maximum value. In particular, we find that stars born with metallicity Z « 0.04 live 
the longest. For larger values of Z, the reduction in nuclear fuel and the change in composition outweigh 
the lifetime extending decrease in luminosity arising from the increased opacity. 

A recent set of galactic chemical evolution calculations (Timmes, 1996) have probed far into the stellif- 
erous epoch. The best indications suggest that the galactic abundance set will approach an asymptotically 
constant composition {X ~ 0.2, Y ~ 0.6, and Z ~ 0.2) over a time scale rj ~ 12. As a consequence, 
any generations of stars formed after 77 ~ 12 will suffer significantly shorter lifetimes than the theoretical 
maximum implied by equation [ 2 . 10 ]. 


2. Stellar Masses vs Metallicity 


The maximum stable stellar mass decreases as metallicity increases. On the main sequence, the 
maximum possible mass is reached when the star’s radiation pressure comes to dominate the thermal (gas) 
pressure within the star. Here, we introduce the usual ansatz that the total pressure at the center of the 
star can be written in the form Pc = Pr + Pg, where the thermal gas pressure is given by the fraction 
Pg = [3Pc and, similarly, Pp = (1 — (3)Pc. Using the ideal gas law for the thermal pressure and the 
equation of state for a gas of photons, we can write the central pressure in the form 


Pc 


[3(l-/3)l 

1/3 

■ kpc ■ 

-a 


-fimp- 


( 2 . 11 ) 


where k is Boltzmann constant and a is the radiation constant. The quantity /i is again the mean molecular 
weight and can be written in the form of equation [2.7]. In hydrostatic equilibrium, the central pressure 
required to support a star of mass M* can be expressed as 


Pc 


-, 1/3 


L36J 


gmV'^ , 


( 2 . 12 ) 
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where pc is the central density (see Phillips, 1994). 

Equating the above two expressions [2.11] and [2.12], we can solve for the mass to find 


M, 


|-108 (1-/3)] 

1/2 

r k 1 

- TTO /3^ - 


.firup. 


G-m 


40Mo , 


(2.13) 


where we have set (3 = 1/2 to obtain the numerical value. The maximum stellar mass thus depends 
somewhat sensitively on the mean molecular weight /x, which in turn is a function of the metallicity. By 
applying the approximations [2.7], [2.8], and [2.9], one can write the maximum mass in the form 

M* = 40M© |(2 - 5Yp/4) - (3 + 5//2)Z/2|^ « 114M© (1 - 2.4^)^ . (2.14) 


Thus, for the expected asymptotic value of the metallicity, Z = 0.2, the maximum mass star is only 
Mh, « 3OM0. 

The continuously increasing metallicity of the interstellar medium will also have implications for low 
mass stars. Higher metallicity leads to more effective cooling, which leads to lower temperatures, which in 
turn favors the formation of less massive stars (e.g., see the recent theory of the initial mass function by 
Adams & Fatuzzo, 1996). The IMF of the future should be skewed even more dramatically in favor of the 
faintest stars. 

The forthcoming metallicity increases may also decrease the mass of the minimum mass main sequence 
star as a result of opacity effects (cf. the reviews of Stevenson, 1991; Burrows & Liebert, 1993). Other 
unexpected effects may also occur. For example, when the metallicity reaches several times the solar value, 
objects with mass M^, = 0.04 M© may quite possibly halt their cooling and contraction and land on the 
main sequence when thick ice clouds form in their atmospheres. Such “frozen stars” would have an effective 
temperature of T* « 273 K, far cooler than the current minimum mass main sequence stars. The luminosity 
of these frugal objects would be more than a thousand times smaller than the dimmest stars of today, with 
commensurate increases in longevity. 


C. The Fate of the Earth and the Sun 


A popular and frequently quoted scenario for the demise of the Earth involves destruction through 
evaporation during the Sun’s asymptotic giant branch (AGB) phase. As the Sun leaves the horizontal 
branch and expands to become an AGB star, its outer radius may swell to such an extent that the 
photospheric radius overtakes the current orbital radius of the Earth. If this state of affairs comes to pass, 
then two important processes will affect the Earth: [1] Evaporation of material due to the extreme heat, 
and [2] Orbital decay through frictional drag. This second process drives the Earth inexorably into the 
giant sun, thereby increasing the efficacy of the evaporation process. Once the earth finds itself inside the 
sun, the time scale for orbital decay is roughly given by the time required for the expiring Earth to sweep 
through its mass. Me, in solar material. This short time interval is given by 


Me 

Pq X^orbit 


50 yr. 


(2.15) 


where p© ~ 10“® g/cm^ is the mass density of solar material at the photosphere. Re ~ 6370 km is the 
radius of the Earth, and Uorbit ~ 30 km/s is the orbital speed. Hence, the demise of the Earth will befall it 
swiftly, even in comparison to the accelerated stellar evolution time scale inherent to the asymptotic giant 
branch. The Earth will be efficiently dragged far inside the sun and vaporized in the fierce heat of the 
stellar plasma, its sole legacy being a small (0.01%) increase in the metallicity of the Sun’s distended outer 
envelope. 

Recent work suggests, however, that this dramatic scene can be avoided. When the sun reaches a 
luminosity of ~ lOOL© on its ascent of the red giant branch, it will experience heavy mass loss through the 



action of strong stellar winds. Mass loss results in an increase in the orbital radii of the planets and can 
help the Earth avoid destruction. However, the actual amount of mass loss remains uncertain; estimates 
are based largely on empirical measurements (see Reimers, 1975), but it seems reasonable that the sun 
will diminish to ~ O.TOMq when it reaches the tip of the red giant branch, and will end its AGB phase as 
a carbon white dwarf with mass ~ 0.5Mq. Detailed stellar evolution calculations for the sun have been 
made by Sackmann, Boothryod, & Kraemer (1993). In their best-guess mass loss scenario, they find that 
the orbital radii for both the Earth and Venus increase sufficiently to avoid being engulfed during the AGB 
phase. Only with a more conservative mass loss assumption, in which the Sun retains 0.83Mq upon arrival 
on the horizontal branch, does the solar radius eventually overtake the Earth’s orbit. 


D. Continued Star Formation in the Galaxy 


Galaxies can only live as long as their stars. Hence it is useful to estimate how long a galaxy can 
sustain normal star formation (see, e.g., Shu, Adams, & Lizano, 1987) before it runs out of raw material. 
One would particularly like to know when the last star forms. 

There have been many studies of the star formation history in both our galaxy as well as other disk 
galaxies (e.g., Roberts, 1963; Larson & Tinsley, 1978; Rana, 1991; Kennicutt, Tamblyn, & Gongdon, 1994; 
hereafter KTG). Although many uncertainties arise in these investigations, the results can be roughly 
summarized as follows. The gas depletion time tr for a disk galaxy is defined to be the current mass in 
gas, Mgas, divided by the star formation rate SFR, i.e.. 


_ Mgas 

~ SFR ■ 


(2.16) 


For typical disk galaxies, this time scale is comparable to the current age of the universe; KTG cite a range 
Tij « 5 - 15 Gyr. The actual time scale for (total) gas depletion will be longer because the star formation 
rate is expected to decrease as the mass in gas decreases. For example, if we assume that the star formation 
rate is proportional to the current mass in gas, we derive a total depletion time of the form 


T = tr ln[Mo/M_F] , 


(2.17) 


where Mq is the initial mass in gas and Mp is the final mass. For typical disk galaxies, the initial gas mass 
is Mo ~ IO^^Mq (see Table 5 of KTG). Thus, if we take the extreme case of Mp = 1 Mq, the total gas 
depletion time is only t « 23t/{ « 120 - 350 Gyr. In terms of cosmological decades, the gas depletion time 
becomes 7725 = 11.1 — 11.5. 

Several effects tend to extend the gas depletion time scale beyond this simple estimate. When stars 
die, they return a fraction of their mass back to the interstellar medium. This gas recycling effect can 
prolong the gas depletion time scale by a factor of 3 or 4 (KTG). Additional gas can be added to the 
galaxy through infall onto the galactic disk, but this effect should be relatively small (cf. the review of 
Rana, 1991); the total mass added to the disk should not increase the time scale by more than a factor of 2. 
Finally, if the star formation rate decreases more quickly with decreasing gas mass than the simple linear 
law used above, then the depletion time scale becomes correspondingly larger. Given these complications, 
we expect the actual gas depletion time will fall in the range 


7715 = 12-14. (2.18) 

Thus, by the cosmological decade 77 « 14, essentially all normal star formation in galaxies will have ceased. 
Goincidentally, low mass M dwarfs have life expectancies that are comparable to this time scale. In other 
words, both star formation and stellar evolution come to an end at approximately the same cosmological 
decade. 

There are some indications that star formation may turn off even more dramatically than outlined 
above. Once the gas density drops below a critical surface density, star formation may turn off completely 
(as in elliptical and SO galaxies). The gas may be heated entirely by its slow accretion onto a central black 
hole. 
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These results indicate that stellar evolution is confined to a reasonably narrow range of cosmological 
decades. It is presumably impossible for stars to form and burn hydrogen before the epoch of recombination 
in the universe (at a redshift 2 ~ 1000 and hence rj ~ 5.5). Thus, significant numbers of stars will exist 
only within the range 

5.5 <77 <14. (2.19) 

The current epoch (77 ~ 10) lies near the center of this range of (logarithmic) time scales. On the other 
hand, if we use a linear time scale, the current epoch lies very near the beginning of the stelliferous era. 


E. The Final Mass Function 


When ordinary star formation and conventional stellar evolution have ceased, all of the remaining 
stellar objects will be in the form of brown dwarfs, white dwarfs, neutron stars, and black holes. One way 
to characterize the stellar content of the universe at this epoch is by the mass distribution of these objects; 
we refer to this distribution as the “Final Mass Function” or FMF. Technically, the Final Mass Function 
is not final in the sense that degenerate objects can also evolve and thereby change their masses, albeit on 
vastly longer time scales. The subsequent evolution of degenerate objects is discussed in detail in §IV. 

Two factors act to determine the FMF: [1] The initial distribution of stellar masses (the initial mass 
function [IMF] for the progenitor stars), and [2] The transformation between initial stellar mass and the 
mass of the final degenerate object. Both of these components can depend on cosmological time. In 
particular, one expects that metallicity effects will tend to shift the IMF toward lower masses as time 
progresses. 


The initial mass function can be specified in terms of a general log-normal form for the mass distribution 
Ip = dN/d\mn, 

1 


In 7 />(ln m) = A — \rn/mc\ | 


( 2 . 20 ) 


where A, me, and (cr) are constants. Throughout this discussion, stellar masses are written in solar units, 
i.e., m = M*/( 1 M 0 ). This general form for the IMF is motivated by the both current theory of star 
formation and by general statistical considerations (Adams & Fatuzzo, 1996; Zinnecker, 1984; Larson, 
1973; Elmegreen & Mathieu, 1983). In addition, this form is (roughly) consistent with observations (Miller 
& Scalo, 1979), which suggest that the shape parameters have the values (cr) « 1.57 and me ~ 0.1 for 
the present day IMF (see also Salpeter, 1955; Scalo, 1986; Rana, 1991). The constant A sets the overall 
normalization of the distribution and is not of interest here. 


For a given initial mass function, we must find the final masses mp of the degenerate objects resulting 
from the progenitor stars with a given mass m. For the brown dwarf range of progenitor masses, m < mn, 
stellar objects do not evolve through nuclear processes and hence mp = m. Here, the scale mp ~ 0.08 is 
the minimum stellar mass required for hydrogen burning to take place. 

Progenitor stars in the mass range mn < m < msN eventually become white dwarfs, where the mass 
scale msN ~ 8 is the minimum stellar mass required for the star to explode in a supernova (note that the 
mass scale msN can depend on the metallicity - see Jura, 1986). Thus, for the white dwarf portion of the 
population, we must specify the transformation between progenitor mass m and white dwarf mass mwD- 
The results of Laughlin et al. (1996) indicate that stars with main sequence masses m < 0.4 will undergo 
negligible mass loss in becoming helium white dwarfs. Unfortunately, this relationship remains somewhat 
ill-defined at higher masses, mostly due to uncertainties in red giant mass loss rates (e.g., see Wood, 1992). 
For the sake of definiteness, we adopt the following transformation between progenitor mass and white 
dwarf mass, 

7Tl 

mwD = —-exp[/?TO], ( 2 . 21 ) 

1 -I- am 

with a = 1.4 and (3 = 1/15. This formula is consistent with the models of Wood (1992) over the appropriate 
mass range and approaches the expected form mwD = m in the low mass limit. 

Stars with large initial masses, m > msN, end their lives in supernova explosions and leave behind a 
neutron star (although black holes can also, in principle, be produced). The mass of the remnant neutron 
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star is expected to be near the Chandrasekhar limit mch ~ 1-4, as confirmed in the case of the binary 
pulsar (Manchester & Taylor, 1977). 

To compute the FMF, one convolves the initial mass function with the transformations from progenitor 
stars to white dwarfs and neutron stars. The Final Mass Function that results is shown in Figure 2. For 
comparison, the initial mass function is also shown (as the dashed curve). Notice that the two distributions 
are similar for masses less than the Chandrasekhar mass (~ I. 4 M 0 ) and completely different for larger 
masses. 

Once the FMF has been determined, one can estimate the number and mass fractions of the various 
FMF constituents. We define Nbd to be the fraction of brown dwarfs by number and Mbd to be the 
fraction of brown dwarfs by mass, with analogous fractions for white dwarfs {Mwd and M.wd) and neutron 
stars {Nns and Mns)- For an IMF of the form [2.20] with present day values for the shape parameters, 
we obtain the following number fractions: 

Afgjy = 0.45, AfwD = 0.55, ATns = 0.0026. (2.22) 

Similarly, for the mass fractions one finds 

Mbd = 0.097, Mwd = 0.88, Mns = 0.024. (2.23) 

Thus, brown dwarfs are expected to be present in substantial numbers, but most of the mass will reside 
in the form of white dwarfs. Neutron stars will make a relatively small contribution to the total stellar 
population. The above values for A/)vs, and Mns were obtained under the assumption that all stars 
m > msN ~ 8 produce neutron stars. In reality, a portion of these high mass stars may collapse to form 
black holes instead, but this complication does not materially affect the basic picture described above. 
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III. DEATH OF THE GALAXY 


We have argued that over the long term, the galaxy will incorporate a large fraction of the available 
baryonic matter into stars. By the cosmological decade ry = 14 — 15, the stellar component of the galaxy 
will be in the form of seemingly inert degenerate remnants. Further galactic activity will involve these 
remnants in phenomena which unfold over time scales ranging from ry ~ 15 — 30. This time period is part 
of what we term the Degenerate Era. 

The course of this long term galactic dynamical evolution is dictated by two generalized competing 
processes. First, in an isolated physical system containing any type of dissipative mechanism (for exam¬ 
ple, gravitational radiation, or extremely close inelastic encounters between individual stars), the system 
must evolve toward a state of lower energy while simultaneously conserving angular momentum. The net 
outgrowth of this process is a configuration in which most of the mass is concentrated in the center and 
most of the angular momentum is carried by small parcels at large radii. (The present day solar system 
presents a good example of this process at work.) Alternatively, a second competing trend occurs when 
collisionless relaxation processes are viable. In a galaxy, distant encounters between individual stars are 
effectively collisionless. Over time, stars tend to be evaporated from the system, the end product of this 
process is a tightly bound agglomeration (perhaps a massive black hole) in the center, containing only a 
fairly small fraction of the total mass. Hence, one must estimate the relative efficiencies of both collisionless 
and dissipative processes in order to predict the final state of the galaxy. This same competition occurs 
for physical systems on both larger scales (e.g., galaxy clusters) and smaller scales (e.g., globular clusters). 

In addition to gravitational radiation and dynamical relaxation, occasional collisions between substellar 
objects - brown dwarfs - provide a channel for continued star formation at a very slow rate. Collisions 
and mergers involving two white dwarfs will lead to an occasional type I supernova, whereas rare impacts 
involving neutron stars will engender even more exotic bursts of energy. Such events are impressive today. 
They will be truly spectacular within the cold and impoverished environment of an evolved galaxy. 

A. Dynamical Relaxation of the Galaxy 

A stellar system such as a galaxy relaxes dynamically because of stellar encounters. The characteristic 
time scale associated with this process in the case of purely stellar systems is well known and can be written 
as 

R N 

T-relax- ^ 12in(^/2). 

where R is the size of the system, v is the typical random velocity, and N is the total number of stars (for 
further discussion, see Lightman & Shapiro, 1978; Shu, 1982; Binney & Tremaine, 1987). The logarithmic 
factor appearing in the denominator takes into account the effects of many small angle deflections of stars 
through distant encounters. The time scale for stars to evaporate out of the system is roughly given by 


'^evap — 100 Trelax ^ 10 yr , 


(3.2) 


where we have used i? = 10 kpc, u = 40 km/s, and N = 10^^ to obtain the numerical result. We thus 
obtain the corresponding estimate 


??evap = 19 -f logio[R/10kpc] -f logio[A^/10 


111 


(3.3) 


Thus, stars escape from the galaxy with a characteristic time scale rj 
Dyson, 1979). 


19 — 20 (see also Islam, 1977; 


The stellar dynamical evolution of the Galaxy is more complicated than the simple picture outlined 
above. First, the galaxy is likely to have an extended halo of dark matter, much of which may be in 
non-baryonic form. Since this dark halo does not fully participate in the dynamical relaxation process, the 
halo tends to stabilize the system and makes the stellar evaporation time scale somewhat longer than the 
simple estimate given above. 
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Other dynamical issues can also be important. In globular clusters, for example, mass segregation 
occurs long before stellar evaporation and binary star heating plays an important (actually dominant) 
role in the long term evolution. On the other hand, equation [3.1] is formally valid only if the stars are 
not bound into binary or triple systems. Binary interaction effects can be important for the long term 
evolution of the stellar component of the galaxy. In particular, the presence of binaries can increase the 
effective interaction cross section and can lead to a variety of additional types of interactions. Both three- 
body interactions and binary-binary interactions are possible. As a general rule, interactions lead to hard 
binaries becoming harder and wide binaries becoming softer or even disrupted (“ionized”). Binaries that 
become sufficiently hard (close) can spiral inwards, become mass transfer systems, and eventually explode 
as supernovae. These effects are just now becoming understood in the context of globular cluster evolution 
(for further discussion of these dynamical issues, see, e.g., Chernoff & Weinberg, 1990; Hut et ah, 1992). 

Galaxies in general, and our galaxy in particular, live in groups or clusters. These larger scale systems 
will also undergo dynamical relaxation processes analogous to those discussed above. However, a more 
immediate issue that can affect our galaxy in the relatively near future is the possibility of merging with 
other galaxies in the local group, in particular Andromeda (M31). The orbits of nearby galaxies have 
been of the subject of much study (e.g., Peebles, 1994), but large uncertainties remain. For the current 
separation between the Milky Way and M31 {d = 0.75 Mpc) and radial velocity {vr = 120 km/s), the two 
galaxies will experience a close encounter at a time At = 6 x 10® yr in the future (i.e., at 77 = 10.2). Whether 
this encounter will lead to a collision/merger or simply a distant passage depends on the tangential velocity 
component, which is not well determined. The models of Peebles (1994) suggest that the distance of closest 
approach will lie in the range 20 - 416 kpc, with more models predicting values near the upper end of this 
range. Thus, more work is necessary to determine whether or not the Milky Way is destined to collide 
with M31 in the relatively near future. 

However, even if our galaxy does not collide with M31 on the first pass, the two galaxies are clearly a 
bound binary pair. The orbits of binary galaxy pairs decay relatively rapidly through dynamical friction 
(e.g., Binney & Tremaine, 1987; Weinberg, 1989). Thus, even if a collision does not occur on the first 
passing, M31 and the Milky Way will not survive very long as individual spiral galaxies. On a time scale 
of approximately 77 = 11 — 12 , the entire local group will coalesce into one large stellar system. 


B. Gravitational Radiation and the Decay of Orbits 

Gravitational radiation acts in the opposite direction: it causes orbits to lose energy and decay so 
that the stars move inward. We first consider the case of a galaxy and its constituent stars. As a given 
star moves through the potential well of a galaxy, its orbit decays through gravitational radiation (e.g., 
Misner, Thorne, & Wheeler, 1973; Weinberg, 1972). The rate of energy loss is proportional to the square 
of the quadrapole moment of the whole system (see also Ohanian & RufRni, 1994). For the case in which 
the galaxy has a large scale quadrapole moment (e.g., a bar), the rate of energy loss from gravitational 
radiation can be written in the simple form 


E \cJ 


(3.4) 


where r = 2'kR/v is the orbit time. For a galaxy, the rotation curve is almost flat with a nearly constant 
velocity v ~ 200 km/s. The time scale tgr for gravitational radiation is thus given by 


tgr = 


2TrR /V 
V V 


-)■ 

c/ 




(3.5) 


where Rq = 10 kpc is a reference length scale for the galaxy. We thus obtain the estimate 


VGR = 24 -k logio[i?/10kpc]. (3.6) 

This time scale corresponds to ~ 10^® orbits around the galactic center. Notice that if the stars are 
radiating incoherently in a nearly smooth potential, the time scale becomes longer by a factor of Mgai /M*, 
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where M* is the mass of the star and Mgai is effective galactic mass. Notice also that gravitational orbital 
decay takes substantially longer than stellar evaporation from the galaxy (see the previous section). Thus, 
the evolution of the galaxy will be dominated by the collisionless process, and hence the majority of stellar 
remnants will be ejected into intergalactic space rather than winding up in the galactic core (see also Islam, 
1977; Dyson, 1979; Rees, 1984). 

Gravitational radiation also causes the orbits of binary stars to lose energy and decay. Of particular 
importance is the decay of binary brown dwarf stars. The eventual coalescence of these systems can lead 
to the formation of a new hydrogen burning star, provided that the mass of the entire system is larger 
than the hydrogen burning limit Mh ~ O.OSMq. The time scale tqd for orbital decay can be written 

2 G3M3 ’ ^ ^ 

where M* is the mass of the the stars and i?o is the initial orbital separation. Inserting numerical values 
and writing the result in terms of cosmological decades, we obtain the result 

T^OD = 19.4 + 41ogio[Ro/(lAU)] - 31ogio[M*/(lM0)]. (3.8) 

This result also applies to planetary orbits (see §III.F below). 

C. Star Formation through Brown Dwarf Collisions 

Once all of the interstellar material has been used up, one viable way to produce additional stars is 
through the collisions of brown dwarfs. These objects have masses too small for ordinary hydrogen burning 
to take place and hence their supply of nuclear fuel will remain essentially untapped. Collisions between 
these substellar objects can produce stellar objects with masses greater than the hydrogen burning limit, 
i.e., stars of low mass. We note that the search for brown dwarfs has been the focus of much observational 
work (see, e.g., Tinney, 1995) and the existence of these objects is now on firm ground (e.g., Golimowski 
et ah, 1995; Oppenheimer et ah, 1995). 


1. Collision Time Scales 


After conventional star formation in the galaxy has ceased, the total number of brown dwarfs in the 
galaxy will be Nq. Although the value of iVo is uncertain and is currently the subject of much current 
research (e.g., see Alcock et ah, 1993; Aubourg et ah, 1993; Tinney, 1995), we expect that Nq is roughly 
comparable to the number of ordinary stars in the galaxy today, Nq ~ 10^^ (see §11.C). The rate T at 
which these brown dwarfs collide is given by 


Nav 1 dN 
^ ~N~dt ’ 


(3.9) 


where N is the number of brown dwarfs in a galaxy with volume V, a is the collision cross section (see 
below), and v is the typical relative velocity. This equation can be integrated to obtain 


N{t) 


Nq 

l + t/rc ’ 


where tc is the characteristic time scale 


TC 


r-i 


1022yr, 


(3.10) 


(3.11) 


or, equivalently. 


rye = 22 + log 


10 


F/(20kpc)= 


- log 


10 


v/ (200km/s) 


(3.12) 
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To obtain this numerical value for the time scale, we have assumed that the collision cross section is given 
by the geometrical cross section of the brown dwarfs; this assumption is justified below. We have also the 
used numerical values V ~ (20kpc)^ and v ~ 200 km/s which are characteristic of the galactic halo. 

The estimate of collision rates given here is somewhat conservative. Nearby stellar encounters can 
lead to the formation of binaries through tidal excitation of modes on the stars (see Press & Teukolsky, 
1977; Lee & Ostriker, 1986). These binaries can eventually decay and thereby lead to additional stellar 
collisions. 

The time scale [3.12] is the time required for the halo population of brown dwarfs to change. Notice 
that this time scale is larger than the evaporation time scale calculated in §III.A. This ordering makes sense 
because distant encounters (which lead to evaporation) must be much more frequent than true collisions. 
For T] < rjc, the collision rate of brown dwarfs for the entire galaxy is given by Ttot = N/tc ~ 10“^^ 
yr“^. The typical outcome of a brown dwarf collision will be the production of a stellar object with mass 
M* ~ O.IM 0 , large enough to burn hydrogen. The stellar (main-sequence) lifetime of such a star is roughly 
2.5 X 10^^ yr. This stellar evolutionary time scale is longer than the time scale on which stars are forming. 
As a result, the galaxy will produce many stars through this process and will contain ~ 100 hydrogen 
burning stars for cosmological decades rj > 14. 

Notice that the time scale for producing stars through brown dwarf collisions is generally much shorter 
than the orbit decay time for brown dwarf binaries. For orbital decay, equation [3.8] implies that rj ~ 22.5 
-|- 41og]^g(i?/lAU). Thus, brown dwarf collisions provide the dominant mechanism for continued star 
formation while the galaxy remains intact. 

2. Collision Cross Sections 

To complete this argument, we must estimate the cross section for colliding brown dwarfs. Consider 
two brown dwarfs with a relative velocity Vrei- For simplicity, we consider the case of equal mass brown 
dwarfs with mass m. The orbital angular momentum of the system is given by 

J = mvreib, (3.13) 

where b is the impact parameter. When the two dwarfs collide and form a composite star of mass ~ 2m, 
the angular momentum can be written 

m = f{2m)R^n, (3.14) 

where R is the stellar radius, is the rotation rate, and / is a numerical constant of order unity which 
depends on the internal structure of the star. We next invoke the constraint that the rotation rate of the 
final state must be less than the break-up speed, i.e.. 


^2 d 2 ^ G{2m) 


< 


R 


(3.15) 


Combining the above results, we obtain a bound on the impact parameter b that can lead to a bound final 
system. We thus obtain 

(3,16) 

^rel 

which can be used to estimate the cross section. 


Using typical numerical values, we find that 6 ~ i? ~ 10^*^ cm, which is roughly comparable to the radius 
of the brown dwarf (e.g., Burrows et ah, 1993). 


3. Numerical Simulations and Other Results 
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In order to illustrate the viability of this collision process, we have done a set of numerical simulations 
using smooth particle hydrodynamics (SPH). We find that collisions between substellar objects can indeed 
form final products with masses greater than the minimum mass required to burn hydrogen. Examples of 
such collisions are shown in Figure 3. In these simulations, density structures from theoretical brown dwarf 
models (Laughlin & Bodenheimer, 1993) are delivered into impact with relative velocity 200 km/s. The 
hydrodynamic evolutionary sequences shown are adiabatic. One expects that the emergent stellar mass 
object will contract toward the main sequence on a Kelvin-Helmholtz time scale and then initiate hydrogen 
burning. 

Finally, we note that white dwarfs will also collide in the galactic halo. As outlined in §11.E, we 
expect roughly comparable numbers of white dwarfs and brown dwarfs at the end of the stelliferous era. 
Although the white dwarfs are actually smaller in radial size, they are more massive and hence have a larger 
gravitational enhancement to their interaction cross section. As a result, the net cross section and hence 
the net interaction rate of white dwarfs should be roughly comparable to that of brown dwarfs (§III.C.l). 
When white dwarfs collide with each other, several different final states are possible, as we discuss below. 

If the two white dwarfs are sufficiently massive, it is possible that the collision product will have a 
final mass which exceeds the Chandrasekhar limit {Mch ~ I. 4 M 0 ) and hence can explode in a supernova. 
Using the final mass function (see §11.E and Figure 2), we estimate that roughly one third of the white 
dwarfs will have masses greater than 0.7Mq and hence only about one tenth of the collisions can possibly 
result in an object exceeding the Chandrasekhar mass limit. The supernova rate from these collisions can 
thus be as large as Fsat ~ 10“^^ yr”^, although it will be somewhat smaller in practice due to inefficiencies. 

The most common type of collision is between two low mass white dwarfs - the final mass function 
peaks at the mass scale M* « 0.13Mq. These low mass objects will have an almost pure helium compo¬ 
sition. If the final product of the collision has a mass larger than the minimum mass required for helium 
burning (Mfje ~ 0.3Mq), then the product star could land on the helium main-sequence (see, e.g., Kip- 
penhahn & Weigert, 1990). In order for the star to burn helium, the collision must be sufficiently energetic 
to impart enough thermal energy into the star; otherwise, the star will become just another helium white 
dwarf. Another possibility exists for collisions between white dwarfs of slightly larger masses. If the prod¬ 
uct of the collision has a mass smaller than the Chandrasekhar mass and larger than the minimum mass 
to burn carbon (O.9M0 < M* < I. 4 M 0 ), the product star could land on the carbon main sequence. Thus, 
this mode of late time star formation can lead to an interesting variety of stellar objects. 


D. The Black Hole Accretion Time 


Large black holes tend to accrete stars and gas and thereby increase their mass. The black hole 
accretion time is the characteristic time scale for a black hole in the center of a galaxy to swallow the rest 
of the galaxy. If we consider collisions of the black hole with stars, and ignore the other processes discussed 
above (gravitational radiation and stellar evaporation), the time for the black hole to absorb the stars in 
the galaxy is given by 


V_ 

av ’ 


(3.18) 


where V = is the volume of the galaxy, v is the typical speed of objects in the galaxy {v ~ 200 km/s), 
and a is the effective cross section of the black hole. As a starting point, we write the cross section in the 
form 


cr = A7ri?|, 


(3.19) 


where A is a dimensionless enhancement factor due to gravitational focusing, and Rs is the Schwarzschild 
radius Rs given by 

Rg = GM/f? . (3.20) 


We thus obtain the time scale 


r = 10 


30 


yr 


M/10®Mo i?/10kpc 


A- 


(3.21a) 


16 



or, equivalently, 

?7accrete = 30 - 2 logig[Af/IO^M©] + 3logig[-R/lOkpc] - logio[A] . (3.21b) 

The time scale ?7accrete ~ 30 is much longer than the time scale for both stellar evaporation and 
gravitational radiation (see also the following section). As a consequence, at these late times, all the stars 
in a galaxy will either have evaporated into intergalactic space or will have fallen into the central black 
hole via gravitational radiation decay of their orbits. Of course, as the black hole mass grows, the accretion 
time scale decreases. Very roughly, we expect ~ 1 — 10% of the stars to fall to the central black hole and 
the remainder to be evaporated; the final mass of the central black hole will thus be Mbh ~ 10® — 10^® 
Mq. 

One can also consider this process on the size scale of superclusters. When rj ~ 30, supercluster-sized 
cosmological density perturbations of length R will have long since evolved to nonlinearity, and will be 
fully decoupled from the Hubble flow. One can imagine an ensemble of ~ 10® — 10^® Mq black holes which 
have descended from dead galaxies and are now roaming freely and hoovering up an occasional remaining 
star in the volume i?®. The characteristic time scale for this process is 


Vaccrete = 33 - 2 logig [M/1O®M0] 3 logig [l?/10Mpc] . (3.22) 


As for the case of the galaxy, however, this straightforward scenario is compromised by additional 
effects. Gravitational radiation will continuously cause the orbits of the black holes to decay, and some 
of them may eventually merge. Stellar encounters with both other stars and with the black holes will 
lead to stellar evaporation from the supercluster sized system. Over the long term, one expects that the 
supercluster will consist of a very large central black hole with a majority of the stars and many of the 
original ~ 10® — 10^® Mq galactic black holes escaping to large distances. In other words, the supercluster¬ 
sized system will behave somewhat analogously to the galaxy, except that it will contain a larger size scale, 
a longer time scale, and two widely disparate mass scales (namely, a stellar mass scale M* ~ IM©, and 
a black hole mass scale Mbh ~ 10® — 10^® Mq). Equipartition effects between the two mass scales will 
come into play, and will drive the galactic black holes toward the center while preferentially ejecting the 
stellar remnants. In principle, this hierarchy can extend up to larger and larger perturbation length scales, 
although the relevant time scales and detailed dynamics become more uncertain as one proceeds with the 
extrapolation. 


E. Annihilation and Capture of Halo Dark Matter 


Galactic halos consist largely of dark matter, much of which may reside in non-baryonic form. Al¬ 
though the nature and composition of this dark matter remains an important open question, one of the 
leading candidates is Weakly Interacting Massive Particles, usually denoted as WIMPs. These particles 
are expected to have masses in the range Mw = 10 — 100 GeV and interact through the weak force and 
gravity only (cf. the reviews of Diehl et ah, 1995; Jungman, Kamionkowski & Griest, 1996; see also the 
recent proposal of Kane & Wells, 1996). Many authors have studied the signatures of WIMP annihilation, 
usually with the hope of finding a detectable signal. One can apply the results of these studies to estimate 
the time scale for the depletion of WIMPs from a galactic halo. 

We first consider the case of direct particle-particle annihilation. Following usual conventions, the rate 
Tw for WIMP annihilation in the halo can be written in the form 


Tw = nwicrv) , 


(3.23) 


where nw is the number density of WIMPs in the halo and (av) is the average value of the annihilation 
cross section times velocity. If WIMPs make up a substantial mass fraction of the galactic halo, their 
number density is expected to be roughly nw ~ I cm“®. The typical velocity of particles in the galactic 
halo is ~200 km/s. Using the most naive dimensional argument, we can estimate the interaction cross 
section as 


a ~ 


5 X 10"®®cm® 


r Mw 1® 
-IGeV. 


(3.24) 
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where Mw is the mass of the particle and Gf is the Fermi constant. The true cross section has additional 
factors which take into account spin dependences, mixing angles, and other model dependent quantities 
(see Diehl et ah, 1995; Jungman et ah, 1996); the form [3.24] is thus highly approximate, but adequate 
for our purposes. We also note that the relic abundance of dark matter particles is determined by the 
interaction cross section; in order for the abundance to be cosmologically significant, the interaction cross 
section must be of order a ~ 10“^^ cm^ (see Kolb & Turner, 1990). 

Putting all of the above results together, we can estimate the time scale tw for the population of 
WIMPs to change, 

TW = p-i =-~ 3 X lO^^yr. (3.25) 

nw\(Jv) 

Thus, in terms of cosmological decades, we obtain the annihilation time scale in the form 

_ DO r 1 r 1 1 r 

Vw 22.5 I'^gio -|pj_ 3 Q 3 1 —3 • (3.26) 

LIO '^^cm'^s Llcm 

It takes a relatively long time for WIMPs to annihilate via direct collisions. In particular, the annihilation 
time scale is much longer than the stellar evaporation time scale (§III.A). 

Another important related effect is the capture of WIMPs by astrophysical objects. The process of 
WIMP capture has been studied for both the Sun (Press & Spergel, 1985; Faulkner & Gilliland, 1985) and 
the Earth (Freese, 1986) as a means of helping to detect the dark matter in the halo (see also Krauss, 
Srednicki, & Wilczek, 1986; Gould, 1987, 1991). Although WIMP capture by the Sun and the Earth 
can be important for dark matter detection, the lifetimes of both (main sequence) stars and planets are 
generally too small for WIMP capture to significantly affect the total population of particles in the galactic 
halo. On the other hand, stellar remnants, in particular white dwarfs, can be sufficiently long lived to have 
important effects. 

In astrophysical objects, WIMPs are captured by scattering off of nuclei. When the scattering event 
leads to a final velocity of the WIMP that is less than the escape speed of the object, then the WIMP has 
been successfully captured. For the case of white dwarfs, we can make the following simple estimate of the 
capture process. The mean free path of a WIMP in matter with white dwarf densities is generally less than 
the radius of the star. In addition, the escape speed from a white dwarf is large, roughly ~ 3000 km/s, 
which is much larger than the velocity dispersion of WIMPs in the halo. As a result, to first approximation, 
most WIMPs that pass through a white dwarf will be captured. The WIMP capture rate F w* by a white 
dwarf is thus given by 

rw* = nwCTWDVrel , (3.27) 

where awD ~ 10 ^® cm^ is the cross sectional area of the white dwarf and Vrei ~ 200 km/s is the relative 
velocity. The capture rate is thus 

Tw* ~ 102®s-i. (3.28) 

With this capture rate, a white dwarf star can consume its weight in WIMPs on a time scale of ~ lO^"* 
yr. The total mass in WIMPs in the halo is expected to be a factor of I-IO times the mass of stars, which 
will be mostly in the form of white dwarfs at these late times (§ 11 .E). As a result, the time scale for white 
dwarfs to deplete the entire halo population of WIMPs via capture is roughly given by 

T ~ 10^®yr or 77 ~ 25. (3.29) 

The actual time scales will depend on the fraction of the galactic halo in non-baryonic form and on the 
properties (e.g., mass) of the particles; these quantities remain unknown at this time. 

The annihilation of halo WIMPs has important consequences for both the galaxy itself and for the 
white dwarfs. Basically, the galaxy as a whole loses mass while the white dwarfs are kept hotter than 
they would be otherwise. The population of captured WIMPs inside the star will build up to a critical 
density at which the WIMP annihilation rate is in equilibrium with the WIMP capture rate (see, e.g., 
Jungman et ah, 1996). Furthermore, most of the annihilation products will be absorbed by the star, 
and the energy is eventually radiated away (ultimately in photons). The net result of this process (along 
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with direct annihilation) is thus to radiate away the mass of the galactic halo on the time scales given 
by equations [3.26] and [3.29]. This process competes with the evaporation of stars through dynamical 
relaxation (§III.A) and the decay of stellar orbits through gravitational radiation (§III.B). 

Since the time scale for WIMP evaporation is much longer than the dynamical time scale, the galaxy 
will adiabatically expand as the halo radiates away. In the outer galaxy, the dark matter in the halo 
dominates the gravitational potential well and hence the stars in the outer galaxy will become unbound 
as the halo mass is radiated away. Since WIMPs do not dominate the potential inside the solar circle, the 
corresponding effects on the inner galaxy are relatively weak. 

The white dwarf stars themselves will be kept hot by this WIMP capture process with a roughly 
constant luminosity given by 

Lwd = Tmw^w* = J~mwnw<^WDVrei ~ 4 x 10~^^Lq , (3.30) 

where T is an efficiency factor (expected to be of order unity) which takes into account the loss of energy 
from the star in the form of neutrinos. With this luminosity, the white dwarf has a surface temperature 
T* « 63 K, where we have assumed a typical white dwarf mass M* = 0.5Mq. As a reference point, we 
note that an entire galaxy of such stars has a total luminosity comparable to that of the sun, Lgai ~ IT©. 
However, most of the radiation will be emitted at infrared wavelengths, A ~ 50/im. 

For completeness, we note that axions provide another viable candidate for the dark matter in the 
galactic halo (see Chapter 10 of Kolb & Turner, 1990). These particles arise from solutions to the strong 
CP problem in quantum chromodynamics (see, e.g., Peccei & Quinn, 1977ab; Weinberg, 1978; Wilczek, 
1978). The coupling of the axion to the photon allows the axion to decay to a pair of photons with a 
lifetime given by 

Ta « 2 X 10^^ yr (wa/leV) (3.31) 

where ruo is the mass of the axion; we have assumed representative values for the remaining particle physics 
parameters. Relic axions with sufficient numbers to contribute to the dark matter budget of the universe 
have masses in the range 10“® eV < ma < 10“^ eV, where the value depends on the production mechanism. 
Using these mass values, we obtain an allowed range of axion decay time scales, 

32 < r/a < 47. (3.32) 


F. The Fate of Planets during Galactic Death 


Planets can be loosely defined as objects that are small enough (in mass) to be supported by ordinary 
Coulomb forces rather than by degeneracy pressure. Over the long term, planets suffer from several 
deleterious processes. They can be vaporized by their evolving parent stars, and their orbits can either 
decay or be disrupted. Barring these more imminent catastrophes, planets will evaporate as their protons 
decay (see §IV.H). 

The theory of general relativity indicates that planetary orbits slowly decay via emission of gravita¬ 
tional radiation (see §III.B). To fix ideas, consider a planet orbiting a star of mass M* at an initial orbital 
radius R. Gravitational radiation drives orbital decay on a time scale given by 


r = 


27ri? 

V \cJ 


2.6 X 10^®yr 



-3 


(3.33) 


or, in terms of cosmological decades. 


= 19.4 + 41ogio[R/lAU] - 31ogio[M,/lMo]. 


(3.34) 


In the interim, planets can be dislodged from their parent stars during encounters and collisions with 
interloping stars. The time scale for these dislocations is given by the time interval required to produce a 
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sufficiently close encounter with another star. Very roughly, if a perturbing star intrudes within a given 
planet’s orbit, then the planet is likely to be entirely ejected from the system. This time scale is given by 


1 

nav ’ 


(3.35) 


where n is the number density of stars (~ 0.1 pc“^ in our galaxy today), v is the relative velocity (~ 100 
km/s), and where the cross section a is determined by the orbital radius of the planet (a « Inserting 

these values, one finds 

T = 1.3 X 

77 = 15.1-21ogio[i?/lAU], (3.37) 

where R is the radius of the planetary orbit. 

Comparing equation [3.33] with equation [3.36], we find that the time scale for gravitational radiation 
is equal to that of stellar encounters for planetary orbits of radius R = 0.2 AU, which is about half the 
radius of the orbit of Mercury in our own solar system. One might guess then, that very close planets, such 
as the recently discovered companion to 51 Pegasus (Mayor & Queloz, 1995; Marcy, Butler, & Williams, 
1996), will eventually merge with their parent stars as a result of radiative orbital decay, while planets with 
larger initial orbits (e.g., the giant planets in our solar system) will be stripped away from their parent 
stars as a consequence of stellar encounters. However, since the time scale for stellar evolution ( 77 * < 14) is 
much shorter than the time scale for orbital decay, close-in planets around solar-type stars will be destroyed 
during the red giant phases long before their orbits are constricted by general relativity. Only the inner 
planets of low mass M dwarfs (which experience no giant phases) will find their fate sealed by gravitational 
radiation. 
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IV. LONG TERM FATE OF DEGENERATE STELLAR OBJECTS 


Brown dwarfs, white dwarfs, neutron stars, and black holes have lifetimes which are not only much 
longer than the current age of the universe {r] = 10 ), but also greatly exceed the expected lifetime of the 
galaxy (77 = 20 — 25). Due to a general lack of urgency, the ultimate fate of these objects has not yet been 
extensively considered. Nevertheless, these objects will not live forever. If the proton is unstable, then 
proton decay will drive the long term evolution of degenerate stellar objects. Black holes are essentially 
unaffected by proton decay, but they gradually dissipate via the emission of Hawking radiation. Both 
proton decay and Hawking radiation yield many interesting astrophysical consequences. In the following 
discussion, we work out the details of these processes (see also Dicus et ah, 1982; Feinberg, 1981). 


A. Proton Decay 


In Grand Unified Theories (GUTs), the proton is unstable and has a finite, albeit quite long, lifetime. 
For example, the proton can decay through the process 


+ 


(4.1) 


and the Feynman diagram for this decay process is shown in Figure 4. Many different additional decay 
channels are possible and the details ultimately depend on the particular theory (e.g., see the reviews of 
Langacker, 1981; Perkins, 1984). In particular, we note that many other decay products are possible, 
including neutrinos. If protons are unstable, then neutrons will also be unstable over a commensurate time 
scale. Free neutrons are of course unstable to /3-decay on a very short time scale (~ 10 minutes); however, 
bound neutrons will be unstable through processes analogous to the proton decay modes (e.g., see Figure 
4). In the present context, the protons and neutrons of interest are bound in “heavy” nuclei (mostly carbon 
and helium) within white dwarfs. 


For the simplest class of GUTs, as illustrated by the decay modes shown in Figure 4, the rate of 
nucleon decay Fp is roughly given by 


Lp 


rrip 


(4.2) 


where mp is the proton mass and 05 is a dimensionless coupling parameter (see, e.g., Langacker, 1981; 
Perkins, 1984; Kane, 1993). The mass scale Mx is the mass of the particle which mediates the baryon 
number violating process. The decay rate should also include an extra numerical factor which takes into 
account the probability that the interacting quarks (which participate in the decay) are in the same place at 
the same time; this numerical factor is less than unity so that the proton lifetime is larger by a corresponding 
factor. To a first approximation, the time scale for proton decay is thus given by 


rp 


lO^^yr 


- Mx 14 

-10i6GevJ ’ 


(4.3) 


where we have taken into account the aforementioned numerical probability factor. The corresponding 
cosmological time scale is 

TIP = 37-b41ogio[Mjf/10i®GeV]. (4.4) 

Notice that this time scale has a very sensitive dependence on the mass scale Mx of the mediating boson. 

We want to find the allowed range for the proton lifetime. This time scale is constrained from below 
by current experimental limits on the lifetime of the proton (e.g., Perkins, 1984). The proton lifetime must 
be greater than 77 ~ 32 (10^^ yr), where the exact limit depends on the particular mode of proton decay 
(Particle Data Group, 1994). Finding an upper bound is more difficult. If we restrict our attention to the 
class of proton decay processes for which equation [4.4] is valid, then we must find an upper bound on the 
mass scale Mx- Following cosmological tradition, we expect the scale Mx to be smaller than the Planck 
scale Mpi « 10^® GeV, which implies the following range for the proton lifetime. 


32 < r7p < 49 . 


(4.5) 
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The lower bound is set by experimental data; the upper bound is more suggestive than definitive (see also 
§IV.F). 

We can find a more restrictive range for the proton lifetime for the special case in which the decay 
mode from some GUT is responsible for baryogenesis in the early universe. (Note that some baryon number 
violating process is necessary for baryogenesis to take place - see Sakharov, 1967). Let us suppose that 
the decay mode from some GUT is valid and that baryogenesis takes place at an energy scale in the 
early universe Eb ~ Mx- This energy scale must be less than the energy scale Ej of the inflationary 
epoch (Guth, 1981). The inflationary energy scale is constrained to be less than ~ 10“^ Mp\ in order to 
avoid overproducing scalar density perturbations and gravitational radiation perturbations (Lyth, 1984; 
Hawking, 1985; Krauss & White, 1992; Adams & Freese, 1995). Gombining these two constraints, we 
obtain the following suggestive range for the time scale for proton decay, 

32 < ryp < 41. (4.6) 

Although a range of nine orders of magnitude in the relevant time scale seems rather severe, the general 
tenor of the following discussion does not depend critically on the exact value. For the sake of definiteness, 
we adopt rjp = 37 as a representative time scale. 

B. White Dwarfs Powered by Proton Decay 

On a sufficiently long time scale, the evolution of a white dwarf is driven by proton decay. When a 
proton decays inside a star, most of the primary decay products (e.g., pions and positrons) quickly interact 
and/or decay themselves to produce photons. For example, the neutral pion 7r° decays into a pair of 
photons with a lifetime of ~ 10“^® sec; positrons, e+, last only ~ 10“^® sec before annihilating with an 
electron and producing gamma rays. Therefore, one common net result of proton decay in a star is the 
eventual production of four photons through the effective reaction 

p + e“^7 + 7 + 7 + 7, (4.7) 

where the typical energy of the photons is given hy E~^ ~ mp/A ~ 235 MeV. These photons have a 
relatively short mean free path within the star and will thermalize and diffuse outwards through a random 
walk process with a characteristic time scale of ~ 10^ yr, much shorter than the evolutionary time scale 
of the system. Additionally, some fraction of the decay products are in the form of neutrinos, which 
immediately leave the system. 

When proton decay is a white dwarf’s primary energy source, the luminosity is 

L*(t) = TNoETpe-'^^^ « TM{t)Tp, (4.8) 

where Nq ~ 10®^ is the initial number of protons in the star, A ~ 1 GeV is the net energy produced per 

decay, and Fp is the decay rate. The factor T is an efficiency parameter which takes into account the 
fraction of energy lost in the form of neutrinos. Very roughly, we expect ~ 1/3 of the energy in the decay 
products to be in neutrinos and hence F ~ 2/3 (e.g., Dicus et ah, 1982). The exact value of the fraction F 
depends on the branching ratios for a particular GUT and hence is model dependent. For a typical decay 
rate of Fp ~ 10“^^ yr”^, the luminosity in solar units becomes 

L* ~ 10-24^0 . (4.9) 

It is perhaps more illuminating to express this stellar luminosity in ordinary terrestrial units. A white 
dwarf fueled by proton decay generates approximately 400 Watts, enough power to run a few light bulbs, or, 
alternately, about 1/2 horsepower. An entire galaxy of such stars has a total luminosity of Lgai ~ 
which is much smaller than that of a single hydrogen burning star. 

The total possible lifetime for a star powered by proton decay is given by 

T = -^ ln[iVo/A^„,i„], (4.10) 
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where Nq ~ 10^^ is the initial number of nucleons in the star and iVmin is the minimum number of nucleons 
required to consider the object a star. If, for example, one takes the extreme case of A^min = !> the time 
required for the star to completely disappear is is < « 130/r p; in general we obtain 


V*=VP + logic [ln(A^o/A^min) • 


(4.11) 


As we show in §1V.D, the object ceases to be a star when A^min ~ 10'^® and hence r]^, r]p + 1.3. 
During the proton decay phase, the stellar surface temperature is given by 


4 ^ TNqETp ppt 

AirapR'i 


(4.12) 


where we have assumed that the spectral energy distribution is simply a blackbody (ap is the Stefan- 
Boltzmann constant). For a IMq star and the typical decay rate Fp, the effective stellar temperature is 
T* ~ 0.06 K. This temperature will be enormously hotter than the temperature of the universe’s background 
radiation at the cosmological decade ry = 37. 

As a white dwarf loses mass via proton decay, the star expands according to the usual mass/radius 
relation 

T 2 

= 0.114- , (4.13) 

GmeiriJ 

where Z and A are the atomic number and atomic weight of the white dwarf material (e.g., Chandrasekhar, 
1939; Shu, 1982; Shapiro & Teukolsky, 1983). For simplicity, we will take typical values and use A = 2Z. 
If we also rewrite the white dwarf mass/radius relation in terms of natural units, we obtain the relation 



(4.14) 


While the white dwarf is in the proton decay phase of its evolution, the star follows a well defined 
track in the H-R diagram, i.e.. 


L* = Lo(T,/To)12/5 , 


(4.15) 


or, in terms of numerical values. 


L* = IQ-^^Lq 


T* 


L0.06K 


12/5 


(4.16) 


We note that the white dwarf mass/radius relation depends on the star’s chemical composition, which 
changes as the nucleons decay (see the following section). This effect will cause the evolutionary tracks to 
depart slightly from the 12/5 power-law derived above. However, this modification is small and will not 
be considered here. 


C. Chemical Evolution in White Dwarfs 

Over the duration of the proton decay phase, the chemical composition of a white dwarf is entirely 
altered. Several different effects contribute to the change in chemical composition. The nucleon decay 
process itself directly alters the types of nuclei in the star and drives the chemical composition toward 
nuclei of increasingly lower atomic numbers. However, pycnonuclear reactions can occur on the relevant 
(long) time scales and build nuclei back up to higher atomic numbers. In addition, spallation interactions 
remove protons and neutrons from nuclei; these free nucleons then interact with other nuclei and lead to 
further changes in composition. 

In the absence of pycnonuclear reactions and spallation, the chemical evolution of a white dwarf is a 
simple cascade toward lower atomic numbers. As protons and neutrons decay, the remaining nuclei become 
correspondingly smaller. Some of the nuclear products are radioactive and will subsequently decay. Given 
the long time scale for proton decay, these radioactive nuclei are extremely short-lived. As a result, only 
the stable isotopes remain. At relatively late times, when the total mass of the star has decreased by a 
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substantial factor (roughly a factor of ten as we show below), almost all of the nuclei left in the star will 
be in the form of hydrogen. 

At high densities and low temperatures, nuclear reactions can still take place, although at a slow rate. 
The quantum mechanical zero point energy of the nuclei allows them to overcome the Coulomb repulsion 
and fuse. In natural units, the nuclear reaction rate can be written in the form 


W = exp 


(4.17) 


where /i is the reduced mass of the nucleus, Rq is the average spacing between nuclei, and a is the fine 
structure constant (see Shapiro & Teukolsky, 1983). A slightly different form for this reaction rate can be 
derived by including anisotropic and electron screening effects (Salpeter & Van Horn, 1969), but the basic 
form is similar. The parameter S{E) is a slowly varying function of energy which takes into account the 
probability of two nuclei interacting given that tunneling has occurred. Specifically, the parameter S is 
related to the cross section a{E) through the relation 


a{E) 


SjE) ^ 


(4.18) 


where T is the tunneling transition probability. The parameter S can be determined either from direct 
experiments or from theoretical calculations (see Shapiro & Teukolsky, 1983; Bahcall, 1989). 

In order to evaluate the time scale for pycnonuclear reactions to occur, one needs to determine the spac¬ 
ing i?o of the nuclei, or, equivalently, the number density of particles. Using the white dwarf mass/radius 
relation, we obtain the result 


IJ.Ro = 2.29A 



4060ATOH■“^/^ 


(4.19) 


where A is average the atomic weight of the nuclei and where we have defined m* = M^,/Mq. 

We can now obtain a rough estimate for the efficiency of pycnonuclear reactions building larger nuclei 
within white dwarfs. As a reference point, we note that for a density of p ~ 10® g cm“®, the time scale 
for hydrogen to fuse into helium is ~ 10® yr (e.g., Shapiro & Teukolsky, 1983; Salpeter & van Horn, 1969), 
which is much shorter than the proton decay time scale. However, the form of equation [4.17] shows that the 
rate of nuclear reactions becomes highly suppressed as the reacting nuclei become larger. The exponential 
suppression factor roughly has the form ~ exp[—where the numerical factor (3 « 22. Thus, as the 
quantity Z increases, the rate of nuclear reactions decreases exponentially. For example, if Z = 6 and 
A = 12 (for carbon), this exponential term is a factor of ~ 10“^®® smaller than that for hydrogen. Because 
of this large exponential suppression, fusion reactions will generally not proceed beyond helium during the 
late time chemical evolution considered here. Thus, the net effect of pycnonuclear reactions is to maintain 
the decaying dwarf with a predominantly helium composition down to a lower mass scale. 

Spallation is another important process that affects the chemical evolution of white dwarf stars during 
the epoch of proton decay. The high energy photons produced through proton decay can interact with 
nuclei in the star. The most common result of such an interaction is the emission of a single free neutron, 
but charged particles (protons), additional neutrons, and gamma rays can also result (e.g., Hubbell, Gimm, 
& Overbo, 1980). The free neutrons will be promptly captured by other nuclei in a type of late time s- 
process (the r-process is of course dramatically irrelevant). The free protons can produce heavier nuclei 
through pycnonuclear reactions, as described above. Both of these mechanisms thus allow heavier elements 
to build up in the star, albeit at a very slow rate and a very low abundance. Thus, the process of spallation 
initially produces free neutrons and protons; but these nucleons are incorporated into other nuclei. As a 
result, the net effect of spallation is to remove nucleons from some nuclei and then give them back to other 
nuclei within the star. The result of this redistribution process is to widen the distribution of the atomic 
numbers (and atomic weights) for the nuclei in the star. 

In order to assess the importance of spallation processes, we must consider the interaction cross section. 
To leading order, the cross section for nuclear absorption of photons is a single “giant resonance” with a 
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peak at about 24 MeV for light nuclei and a width in the range F = 3 - 9 MeV. The relative magnitude 
of this resonance feature is ~ 20 mb (see, e.g., Hubbell, Gimm, & Overbo, 1980; Brune & Schmidt, 1974), 
roughly a factor of 30 smaller than the total interaction cross section (which is dominated by scattering 
and pair production). For each proton decay event, ~ 940 MeV of matter is converted into photons, with 
some neutrino losses. When these photons cascade downward in energy through the resonance regime (at 
~ 24 MeV), there will be 20 - 40 photons and about one in 30 will produce a spallation event. Hence, on 
average, each proton decay event leads to approximately one spallation event. 

Spallation products allow the interesting possibility that a CNO cycle can be set up within the star. 
The time scale for pycnonuclear reactions between protons (produced by spallation) and carbon nuclei is 
short compared to the proton decay time scale. The time scale for pycnonuclear reactions between protons 
and nitrogen nuclei is comparable to the proton decay time scale. Thus, in principle, the white dwarf 
can set up a CNO cycle analogous to that operating in upper-main-sequence stars (see Clayton, 1983; 
Kippenhahn & Weigert, 1990; Shu, 1982). The energy produced by this cycle will be small compared 
to that produced by proton decay and hence this process does not actually affect the luminosity of the 
star. However, this cycle will affect the chemical composition and evolution of the star. As usual, the net 
effect of the CNO cycle is to build four free protons into a helium nucleus and to maintain an equilibrium 
abundance of the intermediate nitrogen and oxygen nuclei. 

In order to obtain some understanding of the chemical evolution of white dwarfs, we have performed a 
simple numerical simulation of the process. Figure 5 shows the results of this calculation for a 1 Mq white 
dwarf with an initial chemical composition of pure carbon ^^C. The simulation assumes that radioactive 
isotopes decay immediately as they are formed through the preferred decay modes. For each proton decay 
event, a spallation event also occurs (see above) and leads to the removal of a nucleon from a random 
nucleus; the spallation products are then assumed to fuse immediately and randomly with other nuclei 
through the s-process and pycnonuclear reactions. The spallation process builds up a small abundance of 
nuclei heavier than the original ^^C, particularly ^^C which has a substantial mass fraction at “early” times. 
The white dwarf evolves through successive phases in which smaller and smaller nuclei are the dominant 
elements by mass fraction. The star never builds up a significant lithium fraction due to the immediate 
fission of newly formed ®Be into ot particles. The star has a broad phase during which '^He dominates the 
composition. When the white dwarf has lost about 60% of its original mass, the hydrogen mass fraction 
begins to predominate. 

D. Final Phases of White Dwarf Evolution 

In the final phases in the life of a white dwarf, the star has lost most of its mass through proton 
decay. When the mass of the star becomes sufficiently small, two important effects emerge: The first effect 
is that degeneracy is lifted and the star ceases to be a white dwarf. The second effect is that the object 
becomes optically thin to its internal radiation produced by proton decay and thus ceases to be a star. In 
the following discussion, we present simple estimates of the mass scales at which these events occur. 

When the star has lost enough of its initial mass to become nondegenerate, most of the nucleons in the 
star will be in the form of hydrogen (see the previous section). A cold star composed of pure hydrogen will 
generally have a thick envelope of molecular hydrogen surrounding a degenerate core of atomic hydrogen. 
As the stellar mass continues to decline through the process of proton decay, the degenerate core becomes 
increasingly smaller and finally disappears altogether. This transition occurs when the degeneracy energy, 
the Coulomb energy, and the self-gravitational energy of the star are all comparable in magnitude; this 
event, in turn, occurs when the central pressure Pc drops below a critical value of roughly a few Megabars 
(Pc ~ 10^^ dyne/cm^). The central pressure in a star can be written in the form 

Pc = f3^, (4.20) 

where /3 is a dimensionless number of order unity. Using the white dwarf mass/radius relation in the form 
of equation [4.13] and setting Z = A = 1, we find the central pressure as a function of stellar mass, 

mi , (4.21) 
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or, equivalently (in cgs units), 


/ M \ 10/3 

Pc~2x 10^^ dyne/cm^ (u^) ' 

Combining these results, we find that the mass scale at which the star becomes nondegenerate is 

given by 

(4.23) 

This mass scale is roughly the mass of a giant planet such as Jupiter (for more detailed discussion of this 
issue, see also Hamada & Salpeter, 1961; Shu, 1982). At this point in its evolution, the star has a radius 
i?* ~ O.li ?0 ~ 7x 10® cm and a mean density of roughly p ~ 1 g/cm^; these properties are also comparable 
to those of Jupiter. As a reference point, notice also that neutral hydrogen atoms packed into a cubic array 
with sides equal to one Bohr radius would give a density of 1.4 g/cvo?. At this transition, a star powered 
by proton decay has luminosity L* « 10“^^ Lq and effective surface temperature T* « 0.0034 K. 

Once the star becomes nondegenerate, it follows new track in the H-R diagram. The expressions for 
the luminosity and surface temperature (see equations [4.8] and [4.12]) remain valid, but the mass/radius 
relation changes. Since the density of matter is determined by Coulomb forces for the small mass scales of 
interest, the density is roughly constant with a value po ~ 1 g/cm^. We can thus use the simple relationship 
M* = 47 rpoi?*/ 3 . Combining these results, we obtain the relation 
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or, in terms of numerical values, 
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(4.25) 


This steep power-law implies that the effective temperature of the star does not change very much during the 
final phases of evolution (the mass has to decrease by 12 orders of magnitude in order for the temperature 
to change by a factor of 10). As a result, effective surface temperatures of order T* ~ 10“^ K characterize 
the final phases of stellar evolution. 


As the star loses mass, it also becomes increasingly optically thin to radiation. As an object becomes 
transparent, it becomes difficult to meaningfully consider the remnant as a star. An object becomes 
optically thin when 

i?*ncr < 1, (4.26) 


where n is the number density of targets and cr is the cross section of interaction between the radiation 
field and the stellar material. In this present context, we must consider whether the star is optically thin 
to both the gamma rays produced by proton decay and also to the internal radiation at longer wavelengths 
characteristic of its bolometric surface temperature. This latter condition is required for the radiation field 
to be thermalized. 


We first consider the conditions for which the star becomes optically thin to the gamma rays (with 
energies ~ 250 MeV) produced by proton decay. Since we are considering the interaction of gamma 
rays with matter, we can write the cross section in the form 


3 mi 


(4.27) 


where (7 is a dimensionless number (of order unity) and (Tt is the Thompson cross section. To a rough 
approximation, the density will be p ~ 1 g/cm^ and hence the number density will have a roughly constant 
value n ~ 10®^ cm“^. Using these values, we find that the “star” will be safely optically thick to gamma 
rays provided its characteristic size is larger than about one meter. In other words, the object must be 
as big as a large rock. These rocks will not, however, look very much like stars. At the extremely low 
bolometric temperatures characteristic of the stellar photospheres at these late times, the wavelength of the 
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photospheric photons will be macroscopic and hence will interact much less strongly than the gamma rays. 
As a result, the spectral energy distribution of these objects will suffer severe departures from blackbody 
spectral shapes. 

In order to consider the optical depth of the star to its internal radiation field, we rewrite the condition 
[4.26] using the relation no = pn, where k is the opacity. As derived above (equation [4.24]), the surface 
temperature is a slowly varying function in this final phase of evolution; as a result, the wavelength of 
photons in the stellar photosphere will be of order A ~ 100 cm. The interaction of this radiation with the 
star depends on the chemical purity and the crystal-grain structure of the stellar material. We can obtain a 
very rough estimate of the opacity by scaling from known astrophysical quantities. For interstellar graphite, 
for example, the opacity at A = 100 /rm is roughly k ~ 1 cm^/g and scales with wavelength according to 
K (X (see Draine & Lee, 1984). We thus estimate that the opacity in the outer layers of the star/rock 
will be K ~ 10“® cm^/g. Thus, in order for the star to be optically thick to its internal radiation, its radius 
must be i?* > 10® cm, which corresponds to a mass scale of 

lO^^g. (4.28) 

All of these values should be regarded as highly approximate. 

^From these results, the ultimate future of white dwarfs, and indeed our own sun, becomes clear: A 
white dwarf emerges from degeneracy as a pure sphere of hydrogen when the mass drops below M* ~ 
10“®Mq. Finally, the remaining object becomes transparent to its own internal radiation when its mass 
dwindles to M* ~ 10^^ g, and at this point it is no longer a star. Stellar evolution thus effectively comes 
to an end. 

Just prior to the conclusion of stellar evolution, the white dwarf experiences about 2000 proton decay 
events per second and hence has a luminosity of L* ~ 1O“®®L0 ~ 4 erg/s, and a temperature T* ~ 10“® 
K. The time at which this transition occurs is given by r ~ 2irp^. 

Given these results, we can now describe the complete evolution of a 1.0 Mq star (e.g., the Sun), 
from its birth to its death. The entire evolution of the such a star in the Hertzsprung-Russell diagram 
is plotted in Figure 6. The star first appears on the stellar birthline (Stabler, 1988) and then follows a 
pre-main sequence track onto the main sequence. After exhausting its available hydrogen, the star follows 
conventional post-main sequence evolution, including red giant, horizontal branch, red supergiant, and 
planetary nebula phases. The star then becomes a white dwarf with mass M* « 0.5Mq and cools along 
a constant radius track. The white dwarf spends many cosmological decades rj = 11 — 25 near the center 
of the diagram (L* = 10^'^ W; T* = 63 K), where the star is powered by annihilation of WIMPs accreted 
from the galactic halo. When the supply of WIMPs is exhausted, the star cools relatively quickly and 
obtains its luminosity from proton decay (L* « 400 W). The star then follows the evolutionary track in 
the lower right part of the diagram (with L* ~ T* ^^^) until mass loss from proton decay causes the star to 
become nondegenerate. The star then becomes a rock-like object supported by Coulomb forces and follows 
a steeper track (with L* ~ T*^) in the H-R diagram until it becomes optically thin. At this point, the 
object ceases to be a star and stellar evolution effectively comes to an end. During its entire lifetime, the 
Sun will span roughly 33 orders of magnitude in luminosity, 9 orders of magnitude in mass, and 8 orders 
of magnitude in surface temperature. 


E. Neutron Stars Powered by Proton Decay 


The evolution of neutron stars powered by proton decay is qualitatively similar to that of white dwarfs. 
Since neutron stars are (roughly) the same mass as white dwarfs, and since proton decay occurs on the 
size scale of an individual nucleon, the luminosity of the neutron star is given by equations [4.8] and [4.9]. 
To leading order, the mass/radius relation for a neutron star is the same as that of white dwarfs with 
the electron mass irie replaced by the neutron mass (see equations [4.13] and [4.14]). Neutron stars are 
thus ~2000 times smaller than white dwarfs of the same mass, and have appropriately warmer surface 
temperatures. Neutron stars undergoing nucleon decay follow a track in the H-R diagram given by 


Lh. = 10-^*Lq 


r n 112/5 
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(4.29) 
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The final phases of the life of a neutron star will differ from the case of a white dwarf. In particular, 
the neutrons in a neutron star come out of degeneracy in a somewhat different manner than the electrons in 
a white dwarf. Within a neutron star, the neutrons exist and do not /3-decay (into protons, electrons, and 
anti-neutrinos) because of the extremely high densities, which are close to nuclear densities in the stellar 
interior. On the exterior, however, every neutron star has a solid crust composed of ordinary matter. As a 
neutron star squanders its mass through nucleon decay, the radius swells and the density decreases. The 
outer layers of the star are less dense than the central regions and hence the outer region will experience 
/3-decay first. Thus, as the mass decreases, neutrons in the outer portion of the star begin to /3-decay into 
their constituent particles and the star must readjust itself accordingly; the net effect is that the crust of 
ordinary matter thickens steadily and moves inwards towards the center. Once the stellar mass decreases 
below a critical value Me*, the crust reaches the center of the star and the transition becomes complete. 
At this point, the star will resemble a white dwarf more than a neutron star. 

This process thus defines a minimum mass neutron star (see Shapiro & Teukolsky, 1983), which is 
roughly characterized by the parameters 

Me* = 0.0925 Mg, PC = 1-55 X 10^4gcm-3, 33* = 164km, (4.30) 

where pc is the central density of the star. It is hard to imagine current-day astrophysical processes which 
produce stellar objects near this limit. The transformation from a neutron star to a white dwarf occurs 
with a time scale given by 

r=-^ln[Mo/Mc*]«^, (4.31) 

ip ip 

where Mq « 1.4Mg is the initial mass of the neutron star. Notice that neutron stars have a possible mass 
range of only a factor of ~ 15, considerably smaller than the mass range available to white dwarfs. 


F. Higher Order Proton Decay 


Not all particle physics theories predict proton decay through the process described above with decay 
rate Tp (equation [4.2] and Figure 4). In theories which do not allow proton decay through this first order 
process, the proton can often decay through second order processes and/or through gravitational effects. 
By a second order process, we mean an interaction involving two protons and/or neutrons, i.e., AB = 2, 
where B is the baryon number. The decay rate for these alternate decay channels is typically much smaller 
than that discussed above. In this section, we discuss the decay rates and time scales for these higher 
order processes (see also Feinberg, Goldhaber, & Steigman, 1978; Wilczek & Zee, 1979; Weinberg, 1980; 
Mohapatra & Marshak, 1980). 


We first consider a class of theories which allow baryon number violation, but do not have the proper 
vertices for direct proton decay {AB = 1). In such theories, proton decay can sometimes take place through 
higher order processes {AB > 1). For example, if the quarks in two nucleons interact as shown in Figure 
7, the decay rate is roughly given by 


Fa 


Q 

4 rup 


(4.32) 


Even for this higher order example, the theory must have the proper vertices for this process to occur. 
We note that some theories forbid this class of decay channels and require AB = 3 reactions in order for 
nucleon decay to take place (e.g., Goity & Sher, 1995; Castano & Martin, 1994). For the example shown 
in Figure 7, the decay rate is suppressed by a factor of {mp/MxY ~ 10®^ relative to the simplest GUT 
decay channel. As a result, the time scale for proton decay through this second order process is roughly 
given by 
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(4.33) 


and the corresponding cosmological time scale is 


= 101 -h 81ogio[Mx/10i®GeV]. 


(4.34) 
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In order for this decay process to take place, the protons involved must be near each other. For the case 
of interest, the protons in white dwarfs are (mostly) in carbon nuclei and hence meet this requirement. 
Similarly, the neutrons in a neutron star are all essentially at nuclear densities. Notice, however, that free 
protons in interstellar or intergalactic space will generally not decay through this channel. 

The proton can also decay through virtual black hole processes in quantum gravity theories (e.g., 
Zel’dovich, 1976; Hawking, Page, & Pope, 1979; Page, 1980; Hawking, 1987). Unfortunately, the time 
scale associated with this process is not very well determined, but it is estimated to lie in the range 

104V<TPijff < IQi^V, (4.35) 

with the corresponding range of cosmological decades 


46 < rjPBH < 169. (4.36) 

Thus, within the (very large) uncertainty, this time scale for proton decay is commensurate with the second 
order GUT processes discussed above. 

We note that many other possible modes of nucleon decay exist. For example, supersymmetric theories 
can give rise to a double neutron decay process of the form shown in Figure 8a (see Goity & Sher, 1995). 
In this case, two neutrons decay into two neutral kaons. Within the context of standard GUTs, decay 
channels involving higher order diagrams can also occur. As another example, the process shown in Figure 
8b involves three intermediate vector bosons and thus leads to a proton lifetime approximately given by 

r;p3 = 165 + 121ogio[Mjf/10i®GeV]. (4.37) 

Other final states are possible (e.g., three pions), although the time scales should be comparable. This 
process (Figure 8b) involves only the most elementary baryon number violating processes, which allow 
interactions of the general form qq qq. As a result, this decay mode is likely to occur even when the 
lower order channels are not allowed. 

Finally, we mention the case of sphalerons, which provide yet another mechanism that can lead to 
baryon number violation and hence proton decay. The vacuum structure of the electroweak theory allows 
for the non-conservation of baryon number; tunneling events between the different vacuum states in the 
theory give rise to a change in baryon number (for further details, see Rajaraman, 1987; Kolb & Turner, 
1990). Because these events require quantum tunneling, the rate for this process is exponentially suppressed 
at zero temperature by the large factor / = exp[47r/Q;w] ~ 10^^^, where aw is the fine structure constant 
for weak interactions. In terms of cosmological decades, the time scale for proton decay through this process 
has the form rjp = rio + 172, where rjo is the natural time scale (for no suppression). Using the light crossing 
time of the proton to determine the natural time scale (i.e., we optimistically take rjo = —31), we obtain 
the crude estimate rjp « 141. Since this time scale is much longer than the current age of the universe, 
this mode of proton decay has not been fully explored. In addition, this process has associated selection 
rules (e.g., ’t Hooft, 1976) that place further limits on the possible events which exhibit nonconservation 
of baryon number. However, this mode of baryon number violation could play a role in the far future of 
the universe. 

To summarize this discussion, we stress that many different mechanisms for baryon number violation 
and proton decay can be realized within modern theories of particle physics. As a result, it seems likely 
that the proton must eventually decay with a lifetime somewhere in the range 

52 <rip< 200, (4.38) 

where the upper bound was obtained by using Mx ~ Mpi in equation [4.37]. 

To put these very long time scales in perspective, we note that the total number Nx of nucleons in 
the observable universe (at the present epoch) is roughly Nx ~ 10^®. Thus, for a decay time of ?7 = 100, 
the expected number Np of nucleons that have decayed within our observable universe during its entire 
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history is far less than unity, Njj ~ 10 The experimental difficulties involved in detecting higher order 
proton decay processes thus become clear. 

If the proton decays with a lifetime corresponding to 77 ~ 100 — 200, the evolution of white dwarfs will 
be qualitatively the same as the scenario outlined above, but with a few differences. Since the evolutionary 
time scale is much longer, pycnonuclear reactions will be much more effective at building the chemical 
composition of the stars back up to nuclei of high atomic number. Thus, stars with a given mass will have 
higher atomic numbers for their constituent nuclei. However, the nuclear reaction rate (equation [4.17]) 
has an exponential sensitivity to the density. As the star loses mass and becomes less dense (according to 
the white dwarf mass/radius relation [4.13, 4.14]), pycnonuclear reactions will shut down rather abruptly. 
If these nuclear reactions stop entirely, the star would quickly become pure hydrogen and proton decay 
through a two body process would be highly suppressed. However, hydrogen tends to form molecules at 
these extremely low temperatures. The pycnonuclear reaction between the two protons in a hydrogen 
molecule proceeds at a fixed rate which is independent of the ambient conditions and has a time scale 
of roughly ry « 60 (see Dyson, 1979, Shapiro & Teukolsky, 1983, and §111.C for simple estimates of pyc¬ 
nonuclear reaction rates). This reaction will thus convert the star into deuterium and helium on a time 
scale significantly shorter than that of higher order proton decay. The resulting larger nuclei can then still 
decay through a second or third order process. We also note that this same mechanism allows for hydrogen 
molecules in intergalactic space to undergo proton decay through a two body process. 


G. Hawking Radiation and the Decay of Black Holes 


Black holes cannot live forever; they evaporate on long time scales through a quantum mechanical 
tunneling process that produces photons and other products (Hawking, 1975). In particular, black holes 
radiate a thermal spectrum of particles with an effective temperature given by 


Tbh 


1 
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(4.39) 


where Mbh is the mass of the black hole. The total life time of the black hole thus becomes 

25607r 

Tbh = - G , (4.40) 

9* 

where determines the total number of effective degrees of freedom in the radiation field. Inserting 
numerical values and scaling to a reference black hole mass of 10 ® Mq, we find the time scale 


Tbh = 10®®yr 


Mbh/IO^Mq , 


(4.41) 


or, equivalently, 


Vbh =83-1- 31ogig[MB///lO®M0]. 


(4.42) 


Thus, even a black hole with a mass comparable to a galaxy {Mbh ~ IO^^Mq) will evaporate through this 
process on the time scale 'qBH ~ 98. One important consequence of this result is that for ry > 100, a large 
fraction of the universe will be in the form of radiation, electrons, positrons, and other decay products. 


H. Proton Decay in Planets 

Planets will also eventually disintegrate through the process of proton decay. Since nuclear reactions 
have a time scale (ry ~ 1500) much longer than that of proton decay and hence are unimportant (see Dyson, 
1979), the chemical evolution of the planet is well described by a simple proton decay cascade scenario (see 
§IV.C). In particular, this cascade will convert a planet initially composed of iron into a hydrogen lattice 
in ~ 6 proton half lives, or equivalently, on a time scale given by 

6 In 2 

T-pianet « “p— « 10®®yr; rypianet « 38. (4.43) 
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This time scale also represents the time at which the planet is effectively destroyed. 

During the epoch of proton decay, planets radiate energy with an effective luminosity given by 

Tplanet = .FMpla„et(f) Tp « Lq f , (4.44) 

L Me 

where Me is the mass of the Earth and where we have used a proton decay lifetime of 10^^ yr. The 
efficiency factor T is expected to be of order unity. Thus, the luminosity corresponds to ~ 0.4 mW. 
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V. LONG TERM EVOLUTION OF THE UNIVERSE 


In spite of the wealth of recent progress in our understanding of cosmology, the future evolution of 
the universe cannot be unambiguously predicted. In particular, the geometry of the universe as a whole 
remains unspecified. The universe can be closed {k = +1; > 1), flat (fc = 0; = 1), or open {k = —1; 

12 < 1). In addition, the contribution of vacuum energy density remains uncertain and can have important 
implications for the long term evolution of the universe. 

A. Future Expansion of a Closed Universe 

If the universe is closed, then the total lifetime of the universe, from Big Bang to Big Crunch, can 
be relatively short in comparison with the characteristic time scales of many of the physical processes 
considered in this paper. For a closed universe with density parameter Oq > 1, the total lifetime tu of the 
universe can be written in the form 


TU = Oo(f2o - , (5.1) 

where Hq is the present value of the Hubble parameter (see, e.g., Peebles, 1993). Notice that, by definition, 
the age r£/ ^ oo as Hq ^ 1. Current cosmological observations suggest that the Hubble constant is 
restricted to lie in the range 50 - 100 km s“^ Mpc“^ (e.g., Riess, Press, & Kirshner, 1995), and hence the 
time scale is restricted to be greater than ~ 10 Gyr. Additional observations (e.g., Loh & Spillar, 
1986) suggest that Oq < 2. Using these results, we thus obtain a lower bound on the total lifetime of the 
universe, 

Tu > 207r Gyr. (5.2) 

In terms of the time variable rj, this limit takes the form 

r]u> 10.8. (5.3) 

This limit is not very strong - if the universe is indeed closed, then there will be insufficient time to allow 
for many of the processes we describe in this paper. 

We also note that a closed universe model can in principle be generalized to give rise to an oscillating 
universe. In this case, the Big Crunch occurring at the end of the universe is really a “Big Bounce” and 
produces a new universe of the next generation. This idea originated with Lemaitre (1933) and has been 
subsequently considered in many different contexts (from Tolman, 1934 to Peebles, 1993). 

B. Density Fluctuations and the Expansion of a Flat or Open Universe 

The universe will either continue expanding forever or will collapse back in on itself, but it is not 
commonly acknowledged that observations are unable to provide a definitive answer to this important 
question. The goal of many present day astronomical observations is to measure the density parameter 
12, which is the ratio of the density of the universe to that required to close the universe. However, 
measurements of 12 do not necessarily determine the long term fate of the universe. 

Suppose, for example, that we can ultimately measure 12 to be some value 12o (either less than or 
greater than unity). This value of 12o means that the density within the current horizon volume has a 
given ratio to the critical density. If we could view the universe (today) on a much larger size scale (we 
can’t because of causality), then the mean density of the universe of that larger size scale need not be the 
same as that which we measure within our horizon today. Let 12big denote the ratio of the density of the 
universe to the critical density on the aforementioned larger size scale. In particular, we could measure 
a value 12o < 1 and have 12big > 1, or, alternately, we could measure 12o > 1 and have 12big < 1. This 
possibility has been discussed at some length by Linde (1988, 1989, 1990). 

To fix ideas, consider the case in which the local value of the density parameter is Uq « 1 and the 
larger scale value is flbig = 2 > 1. (Note that 12 is not constant in time and hence this value refers to the 
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time when the larger scale enters the horizon.) In other words, we live in an apparently flat universe, which 
is actually closed on a larger scale. This state of affairs requires that our currently observable universe lies 
within a large scale density fluctuation of amplitude 

Ap _ ftp - ffbig _ 1 

p f^big 2 

where the minus sign indicates that we live in a locally underdense region. Thus, a density perturbation with 
amplitude of order unity is required; furthermore, as we discuss below, the size scale of the perturbation 
must greatly exceed the current horizon size. 

On size scales comparable to that of our current horizon, density fluctuations are constrained to be 
quite small {Apjp ~ 10“®) because of measurements of temperature fluctuations in the cosmic microwave 
background radiation (Smoot et ah, 1992; Wright et ah, 1992). On smaller size scales, additional measure¬ 
ments indicate that density fluctuations are similarly small in amplitude (e.g., Meyer, Cheng, & Page, 1991; 
Gaier et ah, 1992; Schuster et ah, 1993). The microwave background also constrains density fluctuations on 
scales larger than the horizon (e.g., Grischuk & Zel’dovich, 1978), although the sensitivity of the constraint 
decreases with increasing size scale A according to the relation ~ (Ahor/A)^, where Ahor is the horizon size. 
Given that density fluctuations have amplitudes of roughly ~ 10“® on the size scale of the horizon today, 
the smallest size scale Ai for which fluctuations can be of order unity is estimated to be 

Ai ~ 300Ahor « 10® Mpc. (5.5) 

For a locally flat universe (flo ~ 1), density fluctuations with this size scale will enter the horizon at a time 
« 3 X lO^to « 3 X lO^”^ yr, or, equivalently, at the cosmological decade 

r7i«17.5. (5.6) 

This time scale represents a lower bound on the (final) age of the universe if the present geometry is 
spatially flat. In practice, the newly closed universe will require some additional time to re-collapse (see 
equation [5.1]) and hence the lower bound on the total age becomes approximately p > 18. 

The situation is somewhat different for the case of an open universe with flo < 1- H the universe is 
open, then the expansion velocity will (relatively) quickly approach the speed of light, i.e., the scale factor 
will expand according to i? oc t (for this discussion, we do not include the possibility that flo = 1 — e, 
where e <C 1, i.e., we consider only manifestly open cases). In this limit, the (comoving) particle horizon 
expands logarithmically with time and hence continues to grow. However, the speed of light sphere - the 
distance out to which particles in the universe are receding at the speed of light - approaches a constant 
in comoving coordinates. As a result, density perturbations on very large scales will remain effectively 
“frozen out” and are thus prevented from further growth as long as the universe remains open. Because 
the comoving horizon continues to grow, albeit quite slowly, the possibility remains for the universe to 
become closed at some future time. The logarithmic growth of the horizon implies that the time scale for 
the universe to become closed depends exponentially on the size scale Ai for which density perturbations 
are of order unity. The resulting time scale is quite long [p ^ 100), even compared to the time scales 
considered in this paper. 

To summarize, if the universe currently has a nearly flat spatial geometry, then microwave background 
constraints imply a lower bound on the total age of universe, p > 18. The evolution of the universe at later 
times depends on the spectrum of density perturbations. If large amplitude perturbations (Ap/p > 1) 
enter the horizon at late times, then the universe could end in a big crunch at some time p > pi = 17.5. 
On the other hand, if the very large scale density perturbations have small amplitude {Ap/p <C 1), then 
the universe can continue to expand for much longer time scales. If the universe is currently open, then 
large scale density perturbations are essentially frozen out. 

C. Inflation and the Future of the Universe 

The inflationary universe scenario was originally invented (Guth, 1981) to solve the horizon problem 
and the flatness problem faced by standard Big Bang cosmology (see also Albrecht & Steinhardt, 1982; 
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Linde, 1982). The problem of magnetic monopoles was also a motivation, but will not be discussed here. 
In addition, inflationary models which utilize “slowly rolling” scalar fields can produce density fluctuations 
which later grow into the galaxies, clusters, and super-clusters that we see today (e.g., Bardeen, Steinhardt, 
& Turner, 1983; Starobinsky, 1982; Guth & Pi, 1982; Hawking, 1982). 

During the inflationary epoch, the scale factor of the universe grows superluminally (usually exponen¬ 
tially with time). During this period of rapid expansion, a small causally connected region of the universe 
inflates to become large enough to contain the presently observable universe. As a result, the observed 
homogeneity and isotropy of the universe can be explained, as well as the observed flatness. In order to 
achieve this resolution of the horizon and flatness problems, the scale factor of the universe must inflate 
by a factor of , where the number of e-foldings Nj ~ 60. At the end of this period of rapid expansion, 
the universe must be re-thermalized in order to become radiation dominated and recover the successes of 
standard Big Bang theory. 

Since the conception of inflation, many models have been produced and many treatments of the 
requirements for sufficient inflation have been given (e.g., Steinhardt & Turner, 1984; Kolb & Turner, 1990; 
Linde, 1990). These constraints are generally written in terms of explaining the flatness and causality of 
the universe at the present epoch. However, it is possible, or even quite likely, that inflation will solve the 
horizon and flatness problems far into the future. In this discussion, we find the number Nj of inflationary 
c-foldings required to solve the horizon and flatness problems until a future cosmological decade rj. 

Since the number of e-foldings required to solve the flatness problem is (usually) almost the same as 
that required to solve the horizon problem, it is sufficient to consider only the latter (for further discussion 
of this issue, see, e.g., Kolb & Turner, 1990; Linde, 1990). The condition for sufficient inflation can be 
written in the form 

1^1 . . 

{HR\ ^ {HR)b ’ ^ ^ 

where the left hand side of the inequality refers to the inverse of the product of the Hubble parameter and 
the scale factor evaluated at the future cosmological decade rj and the right hand side refers to the same 
quantity evaluated at the beginning of the inflationary epoch. 

The Hubble parameter at the beginning of inflation takes the form 


_ 87r Mf 


(5.8) 


where Mj is the energy scale at the start of inflation (typically, the energy scale M/ 
corresponds to cosmological decade rji ~ —44.5). Similarly, the Hubble parameter at 
can be written in the form 


h: 


Sir 

Ym|^’ 


~ 10^® GeV, which 
some future time r] 

(5.9) 


where the energy scale M^i is defined by 


Piv) = M^ = PoR^ ^ ■ 


(5.10) 


In the second equality, we have written the energy density in terms of its value po at the present epoch and 
we assume that the universe remains matter dominated. We also assume that the evolution of the universe 
is essentially adiabatic from the end of inflation (scale factor Rend) until the future epoch of interest (scale 
factor Rff), i.e.. 


Rend 



fMi ’ 


(5.11) 


where = Tq/Rjj is the GMB temperature at time rj and Tq « 2.7 K is the GMB temperature today. The 
quantity fMj is the GMB temperature at the end of inflation, after thermalization, and we have introduced 
the dimensionless factor / < 1. 
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Combining all of the above results, we obtain the following constraint for sufficient inflation, 


Nj _ Rend Mj Tq Hrj 


(5.12) 


Next, we write the present day energy density po in terms of the present day CMB temperature Tq, 


Po = 


0 5 


(5.13) 


where /? « 100. The number of e-foldings is thus given by 


Ni = ln[Rend/RB] = HMi/PTo] + i Ini?^ - In / . (5.14a) 

Inserting numerical values and using the definition [1.1] of cosmological decades, we can write this constraint 
in the form 

Ni^ 61+ \n[Mi/{10^^GeV)] +i(r?-10) In 10. (5.14b) 

o 

For example, in order to have enough inflation for the universe to be smooth and flat up to the cosmological 
decade rj = 100, we require Nj « 130 e-foldings of inflation. This value is not unreasonable in that Nj = 
130 is just as natural from the point of view of particle physics as the Nj = 61 value required by standard 
inflation. 

We must also consider the density perturbations produced by inflation. All known models of inflation 
produce density fluctuations and most models predict that the amplitudes are given by 


Ap _ 1 

~ 10 ”$“’ 


(5.15) 


where H is the Hubble parameter and 4) is the scalar field responsible for inflation (Starobinsky, 1982; 
Guth & Pi, 1982; Hawking, 1982; Bardeen, Steinhardt, & Turner, 1983). In models of inflation with more 
than one scalar field (e.g., La & Steinhardt, 1989; Adams & Freese, 1991), the additional fields can also 
produce density fluctuations in accordance with equation [5.15]. 

In order for these density fluctuations to be sufficiently small, as required by measurements of the 
cosmic microwave background, the potential H(4>) for the inflation field must be very flat. This statement 
can be quantified by defining a “fine-tuning parameter” Aft through the relation 


Aft 


AH 

(A4>)4 ’ 


(5.16) 


where AH is the change in the potential during a given portion of the inflationary epoch and A<I> is the 
change in the scalar field over the same period (Adams, Freese, & Guth, 1991). The parameter Aft is 
constrained to less than ~ 10“® for all models of inflation of this class and is typically much smaller. 
Aft ^ 10~i2, for specific models. The required smallness of this parameter places tight constraints on 
models of inflation. 

The aforementioned constraints were derived by demanding that the density fluctuations (equation 
[5.15]) are sufficiently small in amplitude over the size scales of current cosmological interest, i.e., from 
the horizon size (today) down to the size scale of galaxies. These density perturbations are generated 
over Ns ^ 8 e-foldings during the inflationary epoch. However, as discussed in §V.B, large amplitude 
density fluctuations can come across the horizon in the future and effectively close the universe (see also 
Linde, 1988, 1989, 1990). In order for the universe to survive (not become closed) up until some future 
cosmological decade p, density fluctuations must be small in amplitude for all size scales up to the horizon 
size at time rj (within an order of magnitude - see equation [5.1]). As a result, inflation must produce small 
amplitude density fluctuations over many more e-foldings of the inflationary epoch, namely 

«8 + i(r/-10)lnl0, (5.17) 
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where ry is the future cosmological decade of interest. For example, for ry = 100 we would require Ns « 77. 
Although this larger value of Ns places a tighter bound on the fine-tuning parameter Xft, and hence a 
tighter constraint on the inflationary potential, such bounds can be accommodated by inflationary models 
(see Adams, Freese, & Guth, 1991 for further discussion). Loosely speaking, once the potential is flat over 
the usual Ns = 8 e-foldings required for standard inflationary models, it is not that difficult to make it flat 
for Ns = 80. 


D. Background Radiation Fields 

Many of the processes discussed in this paper will produce background radiation fields, which can be 
important components of the universe (see, e.g., Bond, Carr, & Hogan, 1991 for a discussion of present 
day backgrounds). Stars produce radiation fields and low mass stars will continue to shine for several more 
cosmological decades (§11). The net effect of WIMP capture and annihilation in white dwarfs (§III.E) will 
be to convert a substantial portion of the mass energy of galactic halos into radiation. Similarly, the net 
effect of proton decay (§IV) will convert the mass energy of the baryons in the universe into radiation. 
Finally, black holes will evaporate as well, (§IV.H), ultimately converting their rest mass into radiation 
fields. As we show below, each of these radiation fields will dominate the radiation background of the 
universe for a range of cosmological decades, before being successively redshifted to insignificance. 

The overall evolution of a radiation held in an expanding universe can be described by the simple 
differential equation, 

^+4|ftad=^(i), (5.18) 

where Prad is the energy density of the radiation held and S{t) is a source term (see, e.g., Kolb & Turner, 
1990). 

Low mass stars will continue to shine far into the future. The source term for this stellar radiation 
can be written in the form 

A*)^) = — , (5.19) 

tuL 

where and n* are the luminosity and number density of the low mass stars. In the second equality, 
we have introduced the present day mass fraction of low mass stars H*, the nuclear burning efficiency 
e* ~ 0.007, the effective steliar lifetime G, and the present day energy density of the universe po- For this 
example, we have written these expressions for a population of stars with only a single mass; in general, 
one should of course consider a distribution of stellar masses and then integrate over the distribution. As 
a further refinement, one could also include the time dependence of the stellar luminosity L* (see §11). 

For a given geometry of the universe, we find the solution for the background radiation field from low 
mass stars, 

Prad* = e* Hh, p{R) f Y ) (5.20) 

where the dimensionless factor / = 1/2 for an open universe and / = 3/5 for a flat universe. This form is 
valid until the stars burn out at time t = t*. After that time, the radiation field simply redshifts in the 
usual manner, prad* ~ 

For the case of WIMP annihilation in white dwarfs, the source term is given by 

Swit) = = HwPoT? ^F , (5.21) 

where L* and n* are the luminosity and number density of the white dwarfs. In the second equality, we 
have written the source in terms of the energy density in WIMPs, where Q.w is the present day mass 
fraction of WIMPs and F is the effective annihilation rate. The solution for the background radiation held 
from WIMP annihilation can be found. 


Pwrb(l) = / p{R) r t , 


(5.22) 
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where the dimensionless factor / is defined above. This form is valid until the galactic halos begin to run 
out of WIMP dark matter at time t ~ ~ 10^^ yr, or until the galactic halo ejects nearly all of its white 

dwarfs. We note that direct annihilation of dark matter will also contribute to the background radiation 
field of the universe. However, this radiation will be highly nonthermal; the annihilation products will 
include gamma rays with characteristic energy ~ 1 GeV. 

For the case of proton decay, the effective source term for the resulting radiation field can be written 

Sp{t) = poR~^Tpe-^^*, (5.23) 

where is the present day contribution of baryons to the total energy density po, Fp is the proton decay 
rate, and E is an efficiency factor of order unity. For a given geometry of the universe, we obtain the 
solution for the background radiation field from proton decay. 


Pprh{t)=^^Bp{R)FiO, 


(5.24) 


where F{^) is a dimensionless function of the dimensionless time variable ^ = F pt. For an open universe, 


whereas for a flat universe, 


F{0 


l-(l + g)e-^ 


(5.25) 


i^(^)=5-2/3 ^^-2/3 (5 26) 

Jo 

where 7(5/3,^) is the incomplete gamma function (Abramowitz & Stegun, 1972). 

For black hole evaporation, the calculation of the radiation field is more complicated because the 
result depends on the mass distribution of black holes in the universe. For simplicity, we will consider a 
population of black holes with a single mass M and mass fraction VIbh (scaled to the present epoch). The 
source term for black hole evaporation can be written in the form 


SBH{t) = ^BH Po R ^7 - 

oTbh 


1 

1 — t/TBH ’ 


(5.27) 


where tbh is the total lifetime of a black hole of the given mass M (see equation [4.37]). For an open 
universe, we obtain the solution for the background radiation field from black hole evaporation 


Pbhr(0 = ^BH P{R) F{^) , 


(5.28) 


where the dimensionless time variable ^ = I/tbh- For an open universe, the dimensionless function F{^) 
is given by 


whereas for a flat universe, 



(5.29) 


(5.30) 


Each of the four radiation fields discussed here has the same general time dependence. For times short 
compared to the depletion times, the radiation fields have the form 

pit)^nxp{R)Txt, (5.31) 

where fix is the present day abundance of the raw material and Fx is the effective decay rate (notice that 
we have neglected dimensionless factors of order unity). After the sources (stars, WIMPs, protons, black 
holes) have been successively exhausted, the remaining radiation fields simply redshift away, i.e., 

p(i) = P(Wd) (i?/i?e„d)“^ (5.32) 
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where the subscript refers to the end of the time period during which the ambient radiation was produced. 

Due to the gross mismatch in the characteristic time scales, each of the radiation fields will provide 
the dominate contribution to the radiation content of the universe over a given time period. This trend 
is illustrated in Figure 9, which shows the relative contribution of each radiation field as a function of 
cosmological time rj. For purposes of illustration, we have assumed an open universe and the following 
source abundances: low mass stars fl* = 10“^, weakly interacting massive particles 0^ = 0.2, baryons 
^1b = 0.05, and black holes ^Ibh = 0.1. At present, the cosmic microwave background (left over from 
the big bang itself) provides the dominant radiation component. The radiation field from star light will 
dominate the background for the next several cosmological decades. At cosmological decade rj ~ 16, the 
radiation field resulting from WIMP annihilation will overtake the starlight background and become the 
dominant component. At the cosmological decade rj ~ 30, the WIMP annihilation radiation field will have 
redshifted away and the radiation field from proton decay will begin to dominate. At much longer time 
scales, rj ~ 60, the radiation field from black hole evaporation provides the dominant contribution (where 
we have used 10® Mq black holes for this example). 

The discussion thus far has focused on the total energy density prad of the background radiation fields. 
One can also determine the spectrum of the background fields as a function of cosmological time, i.e., one 
could follow the time evolution of the radiation energy density per unit frequency. In general, the spectra 
of the background radiation fields will be non-thermal for two reasons: 

[1] The source terms are not necessarily perfect blackbodies. The stars and black holes themselves pro¬ 
duce nearly thermal spectra, but objects of different masses will radiate like blackbodies of different 
temperatures. One must therefore integrate over the mass distribution of the source population. It 
is interesting that this statement applies to all of the above sources. For the first three sources (low 
mass stars, white dwarfs radiating WIMP annihilation products, and white dwarfs powered by proton 
decay), the mass distribution is not very wide and the resulting composite spectrum is close to that 
of a blackbody. For the case of black holes, the spectrum is potentially much wider, but the mass 
distribution is far more uncertain. 

[2] The expansion of the universe redshifts the radiation field as it is produced and thereby makes the 
resultant spectrum wider than a thermal distribution. However, due to the linear time dependence of 
the emission (equation [5.31]), most of the radiation is emitted in the final cosmological decade of the 
source’s life. The redshift effect is thus not as large as one might naively think. 

To summarize, the radiation fields will experience departures from a purely thermal distribution. However, 
we expect that the departures are not overly severe. 

The above results, taken in conjunction with our current cosmological understanding, imply that it is 
unlikely that the universe will become radiation dominated in the far future. The majority of the energy 
density at the present epoch is (most likely) in the form of non-baryonic dark matter of some kind. A 
substantial fraction of this dark matter resides in galactic halos, and some fraction of these halos can be 
annihilated and hence converted into radiation through the white dwarf capture process outlined in §111.E. 
However, an equal or larger fraction of this dark matter resides outside of galaxies and/or can escape 
destruction through evaporation from galactic halos. Thus, unless the dark matter particles themselves 
decay into radiation, it seems that enough non-baryonic dark matter should survive to keep the universe 
matter dominated at all future epochs; in addition, the leftover electrons and positrons will help prevent 
the universe from becoming radiation dominated (see also Page & McKee, I98Iab). 


E. Possible Eflfects of Vacuum Energy Density 

If the universe contains a nonvanishing contribution of vacuum energy to the total energy density, then 
two interesting long term effects can arise. The universe can enter a second inflationary phase, in which 
the universe expands superluminally (Guth, 1981; see also Albrecht & Steinhardt, 1983; Linde, 1982). 
Alternately, the vacuum can, in principle, be unstable and the universe can tunnel into an entirely new 
state (e.g., Coleman, 1977, 1985). Unfortunately, the contribution of the vacuum to the energy density of 
the universe remains unknown. In fact, the “natural value” of the vacuum energy density appears to be 
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larger than the cosmologically allowed value by many orders of magnitude. This discrepancy is generally 
known as the “cosmological constant problem” and has no currently accepted resolution (see the reviews 
of Weinberg, 1989; Carroll, Press, & Turner, 1992). 


1. Future Inflationary Epochs 

We first consider the possibility of a future inflationary epoch. The evolution equation for the universe 
can be written in the form 

(R\‘^ SttG. ^ /r 

^ (PM + Pvac) ) (5.33) 

where R is the scale factor, pM is the energy density in matter, and pvac is the vacuum energy density. 
We have assumed a spatially flat universe for simplicity. The matter density varies with the scale factor 
according to pM ~ R~^, whereas the vacuum energy density is constant. We can define the ratio 


^ — Pv&cl PO j 


(5.34) 


i.e., the ratio of the vacuum energy density to that of the matter density po at the present epoch. We 
can then integrate equation [5.6] into the future and solve for the time Gac at which the universe becomes 
vacuum dominated. We find the result 


t 


vac 


to + T 


sinh ^[1 ]—sinh 


(5.35) 


where to is the present age of the universe and we have defined t = (67rGpo) both time scales to and 
T are approximately 10^° yr. 

Several results are immediately apparent from equation [5.35]. If the vacuum energy density provides 
any appreciable fraction of the total energy density at the present epoch (in other words, if v is not too 
small), then the universe will enter an inflationary phase in the very near future. Furthermore, almost any 
nonvanishing value of the present day vacuum energy will lead the universe into an inflationary phase on 
the long time scales considered in this paper. For small values of the ratio v, the future inflationary epoch 
occurs at the cosmological decade given by 


^inflate « 10 + ^ logig [^] . (5.36) 

For example, even for a present day vacuum contribution as small as ~ 10“^°, the universe will enter an 
inflationary phase at the cosmological decade r?inflate ~ 30, long before protons begin to decay. In other 
words, the traditional cosmological constant problem becomes even more severe when we consider future 
cosmological decades. 

If the universe enters into a future inflationary epoch, several interesting consequences arise. After a 
transition time comparable to the age of the universe at the epoch [5.36], the scale factor of the universe 
will begin to grow superluminally. Because of this rapid expansion, all of the astrophysical objects in the 
universe become isolated and eventually become out of causal contact. In other words, every given co¬ 
moving observer will see an effectively shrinking horizon (the particle horizon does not actually get smaller, 
but this language has become common in cosmology - see Ellis & Rothman, 1993 for further discussion of 
horizons in this context). In particular, astrophysical objects, such as galaxies and stars, will cross outside 
the speed-of-light sphere and hence disappear from view. For these same astrophysical objects, the velocity 
relative to the observer becomes larger than the speed of light and their emitted photons are redshifted to 
infinity. 

2. Tunneling Processes 

We next consider the possibility that the universe is currently in a false vacuum state. In other words, 
a lower energy vacuum state exists and the universe can someday tunnel to that lower energy state. This 
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problem, the fate of the false vacuum, was first explored quantitatively by Voloshin et al. (1974) and 
by Coleman (1977). Additional effects have been studied subsequently, including gravity (Coleman & De 
Luccia, 1980) and finite temperature effects (e.g., Linde, 1983). 

To obtain quantitative results, we consider an illustrative example in which the vacuum energy density 
of the universe can be described by the dynamics of a single scalar field. Once a field configuration becomes 
trapped in a metastable state (the false vacuum), bubbles of the true vacuum state nucleate in the sea of 
false vacuum and begin growing spherically. The speed of the bubble walls quickly approaches the speed of 
light. The basic problem is to calculate the tunneling rate (the decay probability) from the false vacuum 
state to the true vacuum state, i.e., the bubble nucleation rate V per unit time per unit volume. For 
tunneling of scalar fields at zero temperature (generally called quantum tunneling), the four-dimensional 
Euclidean action S '4 of the theory largely determines this tunneling rate. The decay probability V can be 
written in the form 

V = Ke -^^, (5.37) 

where AT is a determinental factor (see Coleman, 1977, 1985). For purposes of illustration, we assume a 
generic quartic potential of the form 

V($) = + b^'^ + c<P + d. (5.38) 

We can then write the action S '4 in the form 


^4 = |^(2-5)-3 77(<5), (5.39) 

where 6 = 8Xb/a^ and where 77 is a slowly varying function which has a value near unity for most of the 
range of possible quartic potentials (Adams, 1993). The composite shape parameter S varies from 0 to 2 
as the potential V ($) varies from having no barrier height to having nearly degenerate vacua (see Figure 
10 ). 

Even though equations [5.37 - 5.39] describe the tunneling rate, we unfortunately do not know what 
potential (if any) describes our universe and hence it is difficult to obtain a precise numerical estimate 
for this time scale. To get some quantitative feeling for this problem, we consider the following example. 
For the case of no tunneling barrier (i.e., for S '4 = 0), the characteristic decay probability is given by 
Vo ^ K ^ My, where My is the characteristic energy scale for the scalar field. For My = 10^® GeV 
(roughly the GUT scale), Vo ~ 10^^® s“^ cm“^. With this decay rate, the universe within a characteristic 
volume My^ would convert from false vacuum to true vacuum on a time scale of ~ 10“^^ s. Clearly, 
however, the actual decay time scale must be long enough that the universe has not decayed by the present 
epoch. In order to ensure that the universe has survived, we require that no nucleation events have occurred 
within the present horizon volume (~ [3000 Mpc]^) during the current age of the universe (~ 10^° yr). This 
constraint implies that the action S '4 must be sufficiently large in order to suppress nucleation, in particular, 

S 4 > 2311nl0 « 532. (5.40) 

The question then becomes: is this value for S 4 reasonable? For the parameter A, a reasonable range of 
values is 0.1 < A < 1; similarly, for 6, we take the range 0.1 < S < 1.9. Using the form [5.39] for the action 
and setting 77 = 1, we find the approximate range 

0.5<S4<3x10^ (5.41) 

Thus, the value required for the universe to survive to the present epoch (equation [5.40]) can be easily 
realized within this simple model. In the future, however, the universe could tunnel into its false vacuum 
state at virtually any time, as soon as tomorrow, or as late as ry = 10^. If and when this tunneling effect 
occurs, the universe will change its character almost completely. The physical laws of the universe, or at 
least the values of all of the physical constants, would change as the phase transition completes (see Sher, 
1989 and Crone & Sher, 1990 for a discussion of changing laws of physics during a future phase transition). 
The universe, as we know it, would simply cease to exist. 
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Vacuum tunneling of the entire universe is certainly one of the more speculative topics considered in 
this paper. Nevertheless, its inclusion is appropriate since the act of tunneling from a false vacuum into a 
true vacuum would change the nature of the universe more dramatically than just about any other physical 
process. 

It is also possible for the universe to spontaneously create “child universes” through a quantum tun¬ 
neling process roughly analogous to that considered above (e.g., Sato et ah, 1982; Hawking, 1987; Blau, 
Guendelman, & Guth, 1987). In this situation, a bubble of false vacuum energy nucleates in an otherwise 
empty space-time. If this bubble is sufficiently large, it will grow exponentially and will eventually be¬ 
come causally disconnected from the original space-time. In this sense, the newly created bubble becomes 
a separate “child universe”. The newly created universe appears quite different to observers inside and 
outside the bubble. Observers inside the bubble see the local universe in a state of exponential expansion. 
Observers outside the bubble, in the empty space-time background, see the newly created universe as a 
black hole that collapses and becomes causally disconnected. As a result, these child universes will not 
greatly affect the future evolution of our universe because they (relatively) quickly become out of causal 
contact. 

One potentially interesting effect of these child universes is that they can, in principle, receive infor¬ 
mation from our universe. Before the newly created universe grows out of causal contact with our own 
universe, it is connected through a relativistic wormhole, which can provide a conduit for information 
transfer and perhaps even the transfer of matter (see Visser, 1995 for further discussion of wormholes and 
transferability). The implications of this possibility are the subject of current debate (for varying points 
of view, see, e.g., Linde, 1988, 1989; Tipler, 1992; Davies, 1994). 

F. Speculations about Energy and Entropy Production in the Far Future 

Thus far in this paper, we have shown that entropy can be generated (and hence work can be done) 
up to cosmological decades rj ~ 100. For very long time scales rj :§> 100, the future evolution of the 
universe becomes highly uncertain, but the possibility of continued entropy production is very important 
(see §VI.D). Here, we briefly assess some of the possible ways for energy and entropy to be generated in 
the far future. 


1. Continued Formation and Decay of Black Holes 


For the case of a flat spatial geometry for the universe, future density perturbations can provide a 
mechanism to produce entropy. These density perturbations create large structures which can eventually 
collapse to form black holes. The resulting black holes, in turn, evaporate by emitting Hawking radiation 
and thus represent entropy (and energy) sources (e.g., see also Page & McKee, 1981a; Frautschi, 1982). 
Density perturbations of increasingly larger size scale A will enter the horizon as the universe continues to 
expand. The corresponding mass scale Mx of these perturbations is given by 

Ma = Mo(^), (5.42) 

where tx is the time at which the perturbation enters the horizon and Mq « 10^^ Mq is the total mass 
within the present day horizon (at time to)- 


The time tx represents the time at which a given perturbation enters the horizon and begins to grow; 
a large structure (such as a black hole) can only form at some later time after the perturbation becomes 
nonlinear. Suppose that a density perturbation has an initial amplitude when it enters the horizon. In 
the linear regime, the perturbation will grow according to the usual relation 

S = Sx[-J , (5.43) 


where S = Ap/p and t > tx (see Peebles, 1993). Using this growth law, the epoch at which the 
perturbation becomes nonlinear can be written in the form 


?7ni = V\- ^ logic • 


(5.44) 
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For example, if the perturbation has an amplitude 6\ = 10“^, then it becomes nonlinear at time ?7ni = 

+ 6. Since we are interested in very long time scales rj > 100, the difference between the horizon crossing 
time rj\ and the time r/ni of nonlinearity is not overly large. 

One possible result of this process is the production of a large black hole with a mass Mbh ~ Mx. 
The time scale for such a black hole to evaporate through the Hawking process is given by 

tjbh = 101 + 37?a , (5.45) 

where we have combined equations [4.42] and [5.42]. Since rjBH ~ ''InO the universe can form black 

holes faster than they can evaporate. Thus, for the case of a geometrically fiat universe, future density 
perturbations can, in principle, continue to produce black holes of increasingly larger mass. In this case, 
the universe will always have a source of entropy - the Hawking radiation from these black holes. 

We note that these bound perturbations need not necessarily form black holes. The material is (most 
likely) almost entirely non-dissipative and collisionless, and will thus have a tendency to form virialized 
clumps with binding energy per unit mass of order ~ 5<P'. Thus, unless the perturbation spectrum is tilted 
so that (5 is of order unity on these much larger scales, the ensuing dynamics is probably roughly analogous 
to that of a a cluster-mass clump of cold dark matter in our present universe. However, even if the mass 
of the entire perturbation does not form a single large black hole, smaller scale structures can in principle 
form black holes, in analogy to those currently in the centers of present-day galaxies. In addition, it is 
possible that the existing black holes can merge faster than they evaporate through the Hawking process 
(see also §111.D). Thus, the possibility remains for the continued existence of black holes in the universe. 

The process outlined here, the formation of larger and larger black holes, can continue as long as the 
universe remains spatially flat and the density perturbations that enter the horizon are not overly large. 
The inflationary universe scenario provides a mechanism to achieve this state of affairs, at least up to 
some future epoch (see §V.C and in particular equation [5.14]). Thus, the nature of the universe in the 
far future r] ^ 100 may be determined by the physics of the early universe (in particular, inflation) at the 
cosmological decade r] ~ —45. 

Notice that at these very late times, r] ^ 100, the matter entering the horizon will already be “pro¬ 
cessed” by the physical mechanisms described earlier in the paper. Thus, the nucleons will have (most 
likely) already decayed and the matter content of the universe will be mostly electrons, positrons, and 
non-baryonic dark matter particles. Annihilation of both e+-e“ pairs and dark matter will occur simul¬ 
taneously with perturbation growth and hence the final mass of the black hole will be less than Mx- This 
issue must be studied in further depth. 

2. Particle Annihilation in an Open Universe 

If the universe is open, however, then future density perturbations are effectively frozen out (see §V.B) 
and the hierarchy of black holes described above cannot be produced. For an open universe, continued 
energy and entropy production is more difficult to achieve. One process that can continue far into the future, 
albeit at a very low level, is the continued annihilation of particles. Electrons and positrons represent one 
type of particle that can annihilate (see also Page & McKee, 1981ab), but the discussion given below 
applies to a general population of particles. 

Consider a collection of particles with number density n. The time evolution of the particle population 
is governed by the simple differential equation 

^ + 3Hn =—{av)ri^, (5.46) 

at 

where H = R/R is the Hubble parameter and {av) is the appropriate average of interaction cross section 
times the speed (e.g., see Kolb & Turner, 1990). Since we are interested in the case for which the expansion 
rate is much larger than the interaction rate, the particles are very far from thermal equilibrium and we 
can neglect any back reactions that produce particles. For this example, we consider the universe to be 
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open, independent of the activity of this particle population. As a result, we can write R (x t and hence 
H = 1/t. We also take the quantity {av) to be a constant in time (corresponding to s-wave annihilation). 
With these approximations, the differential equation [5.46] can be integrated to obtain the solution 


n{t) = ni 



+ Aoo[l — , 


where we have defined the quantity 

Aoo = -niti{av ), 

and where we have invoked the boundary condition 

n{ti) = ni = constant. 


(5.47) 


(5.48) 


(5.49) 


Analogous solutions for particle annihilation can be found for the case of a flat universe {H = 2/it) and 
an inflating universe {H = constant). 

The difference between the solution [5.47] and the simple adiabatic scaling solution n(t) = ni(ti/t)^ 
is due to particle annihilation, which is extremely small but non-zero. This statement can be quantified 
by defining the fractional difference A between the solution [5.47] and the adiabatic solution, i.e., 

A(t) = ^(<) = Aoo[l-(tiA)2]. (5.50) 

Over the entire (future) lifetime of the universe, the comoving fraction of particles that annihilate is given 
by the quantity Aqo, which is both finite and typically much less than unity. For example, if we consider 
the largest possible values at the present epoch (cr « ctt « 10“^^ cm^, m « 10“® cm“^, « 3 x 10^'^s, 

and V = c), then Aoo ~ 10“^. The fraction Aoo will generally be much smaller than this example. The 
fact that the fraction Aoo is finite implies that the process of particle annihilation can provide only a finite 
amount of energy over the infinite time interval rji < ry < oo. 


3. Formation and Decay of Positronium 

Another related process that will occur on long time scales is the formation and eventual decay of 
positronium. This process has been studied in some detail by Page & McKee (1981ab; see also the discussion 
of Barrow & Tipler, 1986); here we briefly summarize their results. The time scale for the formation of 
positronium in a flat universe is given by 


?7form « 85 -f 2{r]p - 37) - ^ logio[f^e], (5.51) 

where rjp is the proton lifetime (see §IV) and where fig is the mass fraction of after proton decay. For 
a flat or nearly flat universe, most of the electrons and positrons become bound into positronium. In an 
open universe, some positronium formation occurs, but most electrons and positrons remain unattached. 

At the time of formation, the positronium atoms are generally in states of very high quantum number 
(and have radii larger than the current horizon size). The atoms emit a cascade of low energy photons 
until they reach their ground state; once this occurs, the positronium rapidly annihilates. The relevant 
time scale for this decay process is estimated to be 

77decay « Ml + 4(ryp - 37) - ^ logio Pe ]. (5.52) 
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VI. SUMMARY AND DISCUSSION 


Our goal has been to present a plausible and quantitative description of the future of the Universe. 
Table I outlines the most important events in the overall flow of time, as well as the cosmological decades 
at which they occur (see equation [1.1]). In constructing this table, representative values for the (often 
uncertain) parameters have been assumed; the stated time scales must therefore be viewed as approximate. 
Furthermore, as a general rule, both the overall future of the universe, as well as the time line suggested in 
Table I, become more and more uncertain in the face of successively deeper extrapolations into time. Some 
of the effects we have described will compete with one another, and hence not all the relevant physical 
processes can proceed to completion. Almost certainly, parts of our current time line will undergo dramatic 
revision as physical understanding improves. We have been struck by the remarkable natural utility of the 
logarithmic “clock”, r;, in organizing the passage of time. Global processes which can characterize the entire 
universe rarely span more than a few cosmological decades, and the ebb and flow of events is dispersed 
quite evenly across a hundred and fifty orders of magnitude in time, i.e., —50 < rj < 100. 


A. Summary of Results 

Our specific contributions to physical eschatology can be summarized as follows: 

[ 1 ] We have presented new stellar evolution calculations which show the long term behavior of very low 
mass stars (see Figure 1). Stars with very small mass (~ O.IMq) do not experience any red giant 
phases. As they evolve, these stars become steadily brighter and bluer, reaching first a maximum 
luminosity, and second, a maximum temperature, prior to fading away as helium white dwarfs. 

[ 2 ] Both stellar evolution and conventional star formation come to an end at the cosmological decade 
rj ~ 14. This time scale only slightly exceeds the longest evolution time for a low mass star. It also 
corresponds to the time at which the galaxy runs out of raw material (gas) for producing new stars. 
The era of conventional stars in the universe is confined to the range 6 < rj < 14. 

[ 3 ] We have introduced the final mass function (FMF), i.e., the distribution of masses for the degenerate 
stellar objects left over from stellar evolution (see Figure 2). Roughly half of these objects will be 
white dwarfs, with most of the remainder being brown dwarfs. Most of the mass, however, will be in 
the form of white dwarfs (see equations [2.22] and [2.23]). 

[ 4 ] We have explored a new mode of continued star formation through the collisions of substellar objects 
(see Figure 3). Although the time scale for this process is quite long, this mode of star formation will 
be the leading source of new stars for cosmological decades in the range 15 < r/ < 23. 

[ 5 ] We have presented a scenario for the future evolution of the galaxy. The galaxy lives in its present state 
until a time of ?7 ~ 14 when both conventional star formation ceases and the smallest ordinary stars 
leave the main sequence. For times rj > 14, the principle mode of additional star formation is through 
the collisions and mergers of brown dwarfs (substellar objects). The galaxy itself evolves through 
the competing processes of orbital decay of orbits via gravitational radiation and the evaporation of 
stars into the intergalactic medium via stellar encounters. Stellar evaporation is the dominant process 
and most of the stars will leave the system at a time rj ~ 19. Some fraction (we roughly estimate 
~0.01-0.10) of the galaxy is left behind in its central black hole. 

[6] We have considered the annhilation and capture of weakly interacting massive particles (WIMPs) in 
the galactic halo. In the absence of other evolutionary processes, the WIMPs in the halo annihilate 
on the time scale r/ ~ 23. On the other hand, white dwarfs can capture WIMPs and thereby deplete 
the halo on the somewhat longer time scale r] ~ 25. The phenomenon of WIMP capture indicates that 
white dwarf cooling will be arrested rather shortly at a luminosity L* ~ I0 “^^Lq. 

[ 7 ] Depending on the amount of mass loss suffered by the Sun when it becomes a red giant, the Earth 
may be vaporized by the Sun during its asymptotic giant phase of evolution; in this case, the Earth 
will be converted to a small (0.01 %) increase in the solar metallicity. In general, however, planets 
can end their lives in a variety of ways. They can be vaporized by their parent stars, ejected into 
interstellar space through close stellar encounters, merge with their parent stars through gravitational 
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radiation, and can eventually disappear as their protons decay. 

[ 8 ] We have discussed the allowed range for the proton lifetime. A firm lower bound on the lifetime 
arises from current experimental searches. Although no definitive upper limit exists, we can obtain 
a suggestive upper “bound” on the proton lifetime by using decay rates suggested by GUTs and by 
invoking the constraint the mass of the mediating boson, Mx < Mp\ ~ 10^® GeV. We thus obtain the 
following expected range for the proton lifetime 

32 < r;p < 49 + 76(Af - 1), (6.1) 

where the integer N is order of the process, i.e., the number of mediating bosons required for the decay 
to take place. Even for the third order case, we have rjp < 201. Quantum gravity effects also lead to 
proton decay with time scales in the range 46 < ryp < 169. Finally, sphalerons imply rjp ~ 140. 

[ 9 ] We have presented a scenario for the future evolution of sun-like stars (see Figure 6 ). In this case, stars 
evolve into white dwarf configurations as in conventional stellar evolution. On sufficiently long time 
scales, however, proton decay becomes important. For cosmological decades in the range 20 < 77 < 35, 
the mass of the star does not change appreciably, but the luminosity is dominated by the energy 
generated by proton decay. In the following cosmological decades, ry = 35 — 37, mass loss plays a 
large role in determining the stellar structure. The star expands as it loses mass and follows the 
usual mass/radius relation for white dwarfs. The chemical composition changes as well (see Figure 5). 
Proton decay by itself quickly reduces the star to a state of pure hydrogen. However, pycnonuclear 
reactions will be sufficient to maintain substantial amounts of helium (^He and ^He) until the mass 
of the star decreases below ~ O.OIMq. During the proton decay phase of evolution, a white dwarf 
follows a well-defined track in the H-R Diagram given by L* oc T* ' . After the stellar mass decreases 
to M* « 1 O“^M 0 , the star is lifted out of degeneracy and follows a steeper track L^, oc in the H-R 
Diagram. 

[ 10 ] If proton decay does not take place through the first order process assumed above, then white dwarfs 
and other degenerate objects will still evolve, but on a much longer time scale. The relevant physical 
process is likely to be proton decay through higher order effects. The time scales for the destruction 
and decay of degenerate stars obey the ordering 

rjP < 'Hbh < riP2 , ( 6 . 2 ) 

where rjp ~ 37 is the time scale for first order proton decay, tjbh ~ 65 is the time scale for a stellar¬ 
sized black hole to evaporate, and r]P 2 ~ 100 — 200 is the time scale for proton decay through higher 
order processes. 

[ 11 ] In the future, the universe as a whole can evolve in a variety of different possible ways. Future density 
perturbations can come across the horizon and close the universe; this effect would ultimately lead 
(locally) to a big crunch. Alternately, the universe could contain a small amount of vacuum energy 
(a cosmological constant term) and could enter a late time inflationary epoch. Finally, the universe 
could be currently in a false vacuum state and hence kevorking on the brink of instability. In this 
case, when the universe eventually tunnels into the true vacuum state, the laws of physics and hence 
the universe as we know it would change completely. 

[ 12 ] As the cosmic microwave background redshifts away, several different radiation fields will dominate 
the background. In the near term, stellar radiation will overtake the cosmic background. Later on, 
the radiation produced by dark matter annihilation (both direct and in white dwarfs) will provide the 
dominant contribution. This radiation field will be replaced by that arising from proton decay, and 
then, eventually, by the radiation field arising from evaporation of black holes (see Figure 9). 


B. Eras of the Future Universe 

Our current understanding of the universe suggests that we can organize the future into distinct eras, 
somewhat analogous to geological eras: 
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[A] The Radiation Dominated Era. —oo < r] < 4. This era corresponds to the usual time period in which 
most of the energy density of the universe is in the form of radiation. 

[B] The Stelliferous Era. 6 < 77 < 14. Most of the energy generated in the universe arises from nuclear 
processes in conventional stellar evolution. 

[C] The Degenerate Era. 15 < r/ < 37. Most of the (baryonic) mass in the universe is locked up 
in degenerate stellar objects: brown dwarfs, white dwarfs, and neutron stars. Energy is generated 
through proton decay and particle annihilation. 

[D] The Black Hole Era. 38 < ry < 100. After the epoch of proton decay, the only stellar-like objects 
remaining are black holes of widely disparate masses, which are actively evaporating during this era. 

[E] The Dark Era. rj > 100. At this late time, protons have decayed and black holes have evaporated. Only 
the waste products from these processes remain: mostly photons of colossal wavelength, neutrinos, 
electrons, and positrons. The seeming poverty of this distant epoch is perhaps more due to the 
difficulties inherent in extrapolating far enough into the future, rather than an actual dearth of physical 
processes. 

C. Experimental and Theoretical Implications 

Almost by definition, direct experiments that test theoretical predictions of the very long term fate of 
the universe cannot be made in our lifetimes. However, this topic in general and this paper in particular 
have interesting implications for present day experimental and theoretical work. If we want to gain more 
certainty regarding the future of the universe and the astrophysical objects within it, then several issues 
must be resolved. The most important of these are as follows: 

[A] Does the proton decay? What is the lifetime? This issue largely determines the fate stellar objects in 
the universe for time scales longer than ry ~ 35. If the proton is stable to first order decay processes, 
then stellar objects in general and white dwarfs in particular can live in the range of cosmological 
decades ry < 100. If the proton is also stable to second order decay processes, then degenerate stellar 
objects can live for a much longer time. On the other hand, if the proton does decay, a large fraction 
of the universe will be in the form of proton decay products (neutrinos, photons, positrons, etc.) for 
times ry > 35. 

[B] What is the vacuum state of the universe? This issue plays an important role in determining the 
ultimate fate of the universe itself. If the vacuum energy density of the universe is nonzero, then the 
universe might ultimately experience a future epoch of inflation. On the other hand, if the vacuum 
energy density is strictly zero, then future (large) densities perturbations can, in principle, enter our 
horizon and lead (locally) to a closed universe and hence a big crunch. 

[C] What is the nature of the dark matter? Of particular importance is the nature of the dark matter 
that makes up galactic halos. The lifetime of the dark matter particles is also of great interest. 

[D] What fraction of the stars in a galaxy are evaporated out of the system and what fraction are accreted 
by the central black hole (or black holes)? This issue is important because black holes dominate the 
energy and entropy production in the universe in the time range 36 < ry < 100 and the mass of a black 
hole determines its lifetime. 

[E] Does new physics occur at extremely low temperatures? As the universe evolves and continues to 
expand, the relevant temperatures become increasingly small. In the scenario outlined here, photons 
from the cosmic microwave background and other radiation fields, which permeate all of space, can 
redshift indefinitely in accordance with the classical theory of radiation. It seems possible that classical 
theory will break down at some point. For example, in an open universe, the CMB photons will have 
a wavelength longer than the current horizon size (~ 3000 Mpc) at a time ry ~ 40, just after proton 
decay. Some preliminary models for future phase transitions have been proposed (Primack & Sher, 
1980; Suzuki, 1988; Sher, 1989), but this issue calls out for further exploration. 
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D. Entropy and Heat Death 


The concept of the heat death of the universe has troubled many philosophers and scientists since 
the mid-nineteenth century when the second law of thermodynamics was first understood (e.g., Helmholz, 
1854; Clausius, 1865, 1868). Very roughly, classical heat death occurs when the universe as a whole reaches 
thermodynamic equilibrium; in such a state, the entire universe has a constant temperature at all points 
in space and hence no heat engine can operate. Without the ability to do physical work, the universe 
“runs down” and becomes a rather lifeless place. Within the context of modern Big Bang cosmology, 
however, the temperature of the universe is continually changing and the issue shifts substantially; many 
authors have grappled with this problem, from the inception of Big Bang theory (e.g., Eddington, 1931) 
to more recent times (Barrow & Tipler, 1978, 1986; Frautschi, 1982). A continually expanding universe 
never reaches true thermodynamic equilibrium and hence never reaches a constant temperature. Classical 
heat death is thus manifestly avoided. However, the expansion can, in principle, become purely adiabatic 
so that the entropy in a given comoving volume of the universe approaches (or attains) a constant value. 
In this case, the universe can still become a dull and lifeless place with no ability to do physical work. We 
denote this latter possibility as cosmological heat death. 

Long term entropy production in the universe is constrained in fairly general terms for a given class of 
systems (Bekenstein, 1981). For a spatially bounded physical system with effective radius R, the entropy 
S of the system has a well defined maximum value. This upper bound is given by 




2ttRE 

he 


(6.3) 


where E is the total energy of the system. Thus, for a bounded system (with finite size i?), the ratio S/E 
of entropy to energy has a firm upper bound. Furthermore, this bound can be actually attained for black 
holes (see Bekenstein, 1981 for further discussion). 

The results of this paper show that cosmological events continue to produce energy and entropy in the 
universe, at least until the cosmological decade rj ~ 100. As a result, cosmological heat death is postponed 
until after that epoch, i.e., until the Dark Era. After that time, however, it remains possible in principle 
for the universe to become nearly adiabatic and hence dull and lifeless. The energy and entropy generating 
mechanisms available to the universe depend on the mode of long term evolution, as we discuss below. 

If the universe is closed (§V.A) or becomes closed at some future time (§V.B), then the universe will end 
in a big crunch and long term entropy production will not be an issue. For the case in which the universe 
remains nearly flat, density perturbations of larger and larger size scales can enter the horizon, grow to 
nonlinearity, and lead to continued production of energy and entropy through the evaporation of black 
holes (see §V.F.I). These black holes saturate the Bekenstein bound and maximize entropy production. 
Cosmological heat death can thus be avoided as long as the universe remains nearly flat. 

On the other hand, if the universe is open, then density fluctuations become frozen out at some finite 
length scale (§V.B). The energy contained within the horizon thus becomes a finite quantity. However, the 
Bekenstein bound does not directly constrain entropy production in this case because the effective size R 
grows without limit. For an open universe, the question of cosmological heat death thus remains open. For 
a universe experiencing a future inflationary phase (§V.E.l), the situation is similar. Here, the horizon is 
effectively shrinking with time. However, perturbations that have grown to nonlinearity will be decoupled 
from the Hubble flow. The largest nonlinear perturbation will thus define a largest length scale A and 
hence a largest mass scale in the universe; this mass scale once again implies a (finite) maximum possible 
amount of energy available to a local region of space. However, the system is not bounded spatially and 
the questions of entropy production and cosmological heat death again remain open. 

To close this paper, we put forth the point of view that the universe should obey a type of Coper- 
nican Time Principle which applies to considerations of the future. This principle holds that the current 
cosmological epoch (ry = 10) has no special place in time. In other words, interesting things can continue 
to happen at the increasingly low levels of energy and entropy available in the universe of the future. 
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FIGURE CAPTIONS 


Figure 1. The Hertzsprung-Russell diagram for low mass stars for time scales much longer than the current 
age of the universe. The labeled curves show the evolutionary tracks for stars of varying masses, from 0.08 
Mq to O.25M0, as well as the brown dwarf track for a substellar object with mass M* = O.OGM©. The 
inset figure shows the main sequence lifetimes as a function of stellar mass. 


Figure 2. The Final Mass Function (FMF) for stars. Solid curve shows the predicted distribution 
m{dN/dm) for the masses of the degenerate stellar objects (brown dwarfs, white dwarfs, and neutron 
stars) remaining at the cosmological epoch when conventional star formation has ceased. The dashed 
curve shows the mass distribution of the initial progenitor population (the initial mass function). 


Figure 3. Numerical simulation of a collision between two brown dwarfs. The two initial objects have 
masses less than that required for hydrogen burning; the final product of the collision is a true star and 
is capable of sustained hydrogen fusion. The two stars collide with a relative velocity of 200 km/s and an 
impact parameter of ~ 1 stellar radius. The top series of panels shows the collision from a side view; the 
bottom series of panels shows the top view. 


Figure 4. Representative Feynman diagrams for proton decay (top diagram) and neutron decay (bottom 
diagram) shown in terms of the constituent quarks (u, d, d). These processes are the form expected for 
the simplest Grand Unified Theories. The particles X and Y are the intermediate vector bosons which 
mediate the baryon number violating process and are expected to have masses comparable to the GUT 
scale ~ 10^® GeV. 


Figure 5. Chemical evolution of a white dwarf star during proton decay. The curves show the mass fractions 
of the major component nuclei in the star as a function of time, which is measured here in terms of the 
stellar mass. The initial state is a 1.0 Mq white dwarf made of pure This simulation includes the 
effects of spallation and radioactivity (see text). 


Figure 6. The the complete evolution of the Sun (or any IMq star) in the H-R Diagram. The track shows 
the overall evolution of a star, from birth to final death. The star first appears in the H-R diagram on the 
stellar birthline and then follows a pre-main sequence track onto the main sequence. After its post-main 
sequence evolution (red giant, horizontal branch, red supergiant, and planetary nebula phases), the star 
becomes a white dwarf and cools along a constant radius track. The star spends many cosmological decades 
ry = 11 — 25 at a point near the center of the diagram (L* = lO^"* W; T* = 63 K), where the star is powered 
by annihilation of WIMPs accreted from the galactic halo. When the supply of WIMPs is exhausted, the 
star cools relatively quickly and obtains its luminosity from proton decay (L* « 400 W). The star then 
follows the evolutionary track in the lower right part of the diagram (with L* ~ until mass loss 

from proton decay causes the star to become optically thin. At this point, the object ceases to be a star 
and stellar evolution comes to an end. 


Figure 7. Representative Feynman diagram for nucleon decay for a AB = 2 process, i.e., a decay involving 
two nucleons. The net result of this interaction (shown here in terms of the constituent quarks) is the 
decay of a neutron and a proton into two pions, n + p —> + Tr~^. The Y particle mediates the baryon 

number violating process. Similar diagrams for neutron-neutron decay and for proton-proton decay can be 
obtained by changing the type of spectator quarks. 

Figure 8. Representative Feynman diagram for higher order nucleon decay processes, shown here in terms 
of the constituent quarks, (a) Double neutron decay for a supersymmetric theory. The net reaction converts 
two neutrons n into two neutral kaons K^. The tildes denote the supersymmetric partners of the particles. 
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(b) Double nucleon decay involving three intermediate vector bosons Y. Other final states are possible 
(e.g., three pions), but the overall decay rate is comparable and implies a decay time scale rjp ~ 165 + 
121ogio[Mv/10i6GeV]. 


Figure 9. Background radiation fields in the universe. The vertical axis represents the ratio of the energy 
density in radiation to the total energy density (assuming the universe remains matter dominated). The 
horizontal axis is given in terms of cosmological decades rj. The various curves represent the radiation 
fields from the cosmic microwave background (CMB), light from low mass stars (S), radiation from WIMP 
annihilation in white dwarfs (WIMPs), radiation from proton decay (p decay), and black hole evaporation 
(black holes). 


Figure 10. Potential G($) of a scalar field which determines the vacuum state of the universe. This 
potential has both a false vacuum state (labeled F) and a true vacuum state (labeled T). As illustrated 
by the dashed curve, the universe can tunnel from the false vacuum state into the true vacuum state at 
some future time. 
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Table I: Important Events in the History and Future of the Universe 


The Big Bang . rj = —oo 

Planck Epoch.-50.5 

GUT Epoch .-44.5 

Electroweak Phase Transition .-17.5 

Quarks become confined into Hadrons.-12.5 

Nucleosynthesis.-6 

Matter Domination .4 

Recombination.5.5 

First possible Stellar Generation .6 

Formation of the Galaxy .9 

Formation of the Solar System .9.5 

Today: The Present Epoch .10 

Our Sun dies .10.2 

Close Encounter of Milky Way with Andromeda (M31) .10.2 

Lower Bound on the Age of closed Universe .10.8 

Lifetime of Main Sequence Stars with Lowest Mass .13 

End of conventional Star Formation.14 

Planets become detached from Stars .15 

Star Formation via Brown Dwarf Collisions .16 

Lower Bound on Age of flat Universe (with future Ap/p >1) .18 

Stars evaporate from the Galaxy .19 

Planetary Orbits decay via Gravitational Radiation .20 

WIMPs in the Galactic Halo annihilate .22.5 

Star Formation via Orbital Decay of Brown Dwarf Binaries.23 

Stellar Orbits in the Galaxy decay via Gravitational Radiation .24 

White Dwarfs deplete WIMPs from the Galactic Halo .25 

Black Holes accrete Stars on Galactic Size Scale .30 

Black Holes accrete Stars on Cluster Size Scale.33 

Protons decay .37 

Neutron Stars /3-decay .38 

Planets destroyed by Proton Decay .38 

White Dwarfs destroyed by Proton Decay .39 
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Axions decay into Photons. 

Hydrogen Molecules experience Pycnonuclear Reactions .60 

Stellar-sized Black Holes evaporate.65 

Black Holes with M = IO^Mq evaporate .83 

Positronium formation in a Flat Universe .85 

Galaxy-sized Black Holes evaporate .98 

Black Hole with Mass of current Horizon Scale evaporates.131 

Positronium decay in a Flat Universe .141 

Higher order Proton Decay Processes .~100 - 200 
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